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ABSTRACT 

This thesis investigates the feasibility of adapting Naval 
Architecture and Marine Engineering design processes to digital computer 
solution. 

Two problems, the determination of critical whirling frequencies 
of propeller shafts and the determination of stresses in a transverse web 
frame are programmed for a digital computer. A detailed explanation of 
the programs and a discussion of the methods of approach and difficulties 
encountered in adapting the solutions of these problems to a digital 
computer is included. The most important results obtained from programming 
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these problems are experience with and knowledge of digital computers. 
Applying this knowledge and experience to the general field of naval 
architecture the following conclusions and recommendations are presented: 



Conclusions : 

a. Digital computers are applicable to naval architecture problems. 

b. The average naval architect has sufficient background to pro- 
gram a digital computer. 

c. A programmer should have a complete engineering understanding 
of the problem and some familiarity with numerical methods of 
mathematical processes. 

d. In general the limitations imposed by the computer on a 
problem are no more restrictive than those imposed ty manual 
methods. 

e. The digital computer used should be able to handle numbers 
usually encountered in engineering problems. 



Recommendations : 

a. Digital computer should be used in naval architecture, 

b. Digital computers are suited for theses and research problems. 

c. Any future thesis should only consider the computer as an aid 
in the solution of a specific problem and not general appli- 
cability of the computer. 

d. For engineering applications the programmer should have; (a) a 
good engineering tinders tanding of his field; (b) knowledge of 
numerical methods of mathematical processes; (c) organizational 
abilities. 



Thesis Supervisor: Robert J. Hansen 

Title: Associate Professor of Structural Engineering 
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Introduction 

The field of naval architecture contains some of the most 
complex problems facing the engineering profession in this modern day 
and age. A ship is a structure containing a network of mechanical 
equipage that is unequaled by any other movable structure in the world. 
The complexity of such a structure almost defies analytic approach and 
as a result much of the design work is based on previously built vessels. 
A parent design is altered by repeating various numerical procedures 
until the desired result is obtained. The practice of naval architecture 
can best be described as laborious and repetitive with a necessity for 
a background of experience. 

In an effort to eliminate the necessity of utilizing personnel 
to perform laborious calculations, following World War II a committee was 
appointed by the Bureau of Ships to investigate and report on the appli- 
cability of utilizing machine aided calculations for problems arising in 
the preliminary design stage. Due to the stage of development of elec- 
tronic computers at that time an analog computer was developed for 
stability calculations in the preliminary design. Although this investi- 
gation did consider the application of digital computers to problems in 
the preliminary design stage, the application to the general field of 
naval architecture was not fully considered, (ll)^ 

The purpose of this thesis is to investigate the applicability 
of digital computers to the field of naval architecture, with primary 
emphasis on structural design problems, to acquaint the authors and other 

1. Numbers in parentheses refer to the Bibliography at the end of the 
thesis. 
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naval architects with this new method of approach, to determine the type 
problems best suited for computer solution and publish some of the diffi- 
culties encountered, to determine the type of person required to utilize 
most effectively the machine aids and to incite interest in applying 
digital computers in the investigation of problems heretofore considered 
too laborious and time consuming to analyze by hand calculations. To 
do this the authors had to first acquaint themselves with an existing 
computer and secondly, to pick sample engineering problems of generic 
nature for specific programming. 

The authors are graduate students in naval architecture whose 
studies have been primarily channeled in the direction of hull structures. 

The background in mathematics included calculus through differential 
equations and a basic course in vector and complex mathematics. Their 
background in mathematics is no better than can be expected from the 
practicing naval architect. 

Two problems were picked as samples, one, calculation of the 
whirling frequencies of a propeller shafting arrangement and two, an 
analysis of an indeterminate structure of five degrees of indeterminancy. 

These problems were picked because they are examples of types that require 
reiterative procedures in their solutions. Further, both problems had 
solutions that had been previously determined by hand. 

The digital computer investigated and utilized was the Whirlwind I 
computer at MIT. This computer started operations in November 1950 and 
has been used as a guide for the development of many commercially available 
machines today. For this reason the general conclusions drawn from an 
analysis of the applicability of Whirlwind I should apply to any machine 
considered for particular problems, although other machines may have different 
capabilities. 



Basically", all digital computers utilize the same theory of 
operation. They all contain the following elements: storage element, 

control element, arithmetic element and input and output devices, as 
depicted in Figure I. 

Calculations are performed in the arithmetic element utilizing 
the numbers and instructions contained in the storage element. Input 
and output devices are used to get information in and out of the computer. 
The control element makes the components function. 




Schematic Diagram of Digital Computer 
Figure I 



The discussion of the individual elements that follows pertains 
to Whirlwind I only, although other computers have similar equipment that 
performs the same functions. 

The storage element consists of 20U8 registers or "cells." 
Instructions or numbers are stored in these registers by impressing 
electrical pulses on magnetic cores. Each register contains 16 such core 
and, therefore, the basic word-length is 16 digits. The pulses may be 
either positive or negative and the resulting direction of magnetization 
determines whether the digit represented is 0 or 1. This means that 
all the machine claculations are carried out in the binary number system. 
In every register the first binary digit represents the sign of the number 
and the remaining 15 binary digits can represent approximately 5 decimal 
digits. 

Instructions are stored in the registers in a slightly different 
manner. Here the first five binary digits represent the mathematical 
operation that is to be performed when the instruction is obeyed. The 
remaining eleven binary digits represent the location of the register 
containing the number on which the operation is to be performed. 

The permanent storage of Whirlwind I may be increased to 
approximately 22,000 registers by the use of a magnetic drum. This device 
is somewhat slower in the transfer of information than is the permanent 
storage, but it is very convenient for the solution of long problems. 

The arithmetic element consists of three registers. In one of 
these, called the accumulator, or "AC" register, the basic arithmetic 
operations of addition and subtraction are performed. Multiplication and 
division are accomplished by a series of these two basic operations. 
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The input devices are a photoelectric tape reader and a mechanical 
tape reader. These devices read into the storage element coded information 
from a punched paper tape. They convert holes in the tape into electrical 
signals to properly magnetize the registers of the storage element. A 
conversion routine is available to transform decimal information on the 
tape into binary information for storage in the computer. This feature 
is a great aid in programming and enhances the flexibility of the computer. 

The output devices available are flexowriters, five magnetic 
tape units and an oscilloscope and camera. The desired output information 
can be typed out directly by flexowriters, or it can be stored on magnetic 
tape and typed out by the flexowriter at a later time. The latter procedure 
saves computer time since it is much faster to read information onto 
magnetic tape than to type it. Another method of getting the information 
out is to display it on the oscilloscope and photograph it. This saves 
typing time, but the film must be processed. However, it is convenient 
for reproducing the results. Standard output programs are available or 
special ones may be written to get the information out in the form desired. 
For instance, numbers may be in octal or decimal system; written in columns, 
lines, blocks, or tabulated; formed raised to a fixed power of ten or with 
a fixed decimal point. This flexibility of output is an important feature 
of digital computers. 



Control Element 

Each of the various components of the computer consists of 
combinations of electrical circuitry and transfer buses. Only through 
intelligent control of these devices can electrical signals be made to 
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perform arithmetic operations. The control element performs this function 
by causing the flow of magnetic signals to reach the proper element at the 
appropriate time. Basically then, the control element is a switching 
device that channels the information through transfer buses to itself 
and the other elements. To do this, the control element must interpret 
the program instructions, which are themselves magnetic signals, into 
other magnetic signals which control the operation of the computer. 

The speed of operation of the computer depends upon the memory 
access time and the time required to integrate a number of switching com- 
binations into a useful operation. Switching is performed by utilizing 
electrical impulses from a high frequency oscillator. A new set of 
switching combinations can be established for each time pulse received 
by the control element. The time required to perform an operation depends 
on the number of impulses necessary to perform the switching and the 
repetition rate of the impulses. In Whirlwind I the pulse repetition rate 
is 1 micro-second. At this rate the addition of two numbers requires 
micro-seconds and multiplication requires approximately 65 micro-seconds. 
Using the conventional rating system whiich is the number of additions a 
computer can perform per second, the Whirlwind I is rated at 20,000 opera- 
tions per second. 

The instructions necessary to perform the desired mathematical 
calculations are called a program and are fully explained in Appendix A. 
For each computer there exists a standard instruction code. By ingenious 
combinations of these standard instructions the programmer can cause a 
computer to perform the requisite numerical calculations. 
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One of the more important instructions is the "conditional 
response." This instruction allows the computer to make an either-or 
decision based on the sirn of the number in the accumulator. As a 
result of this capability the computer can be made to cycle back through 
a series of instructions any desired number of times. The cycling ability 
is ideally suited to repetitive processes which are useful in design calcu- 
lations. 

Digital computers have been widely used in other fields of 
engineering and new adaptations are being developed every day. Some of 
the more important uses developed are: 

(1) Aircraft weight, balance and thrust calculations. 

(2) Steam turbine design 

(3) Catapault launching analysis 

(li) General structural application in the aircraft 
industry. 
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PROCEDURE 



General 



In order to accomplish the objectives discussed in the 
introduction it was necessary that some knowledge of computers be obtained. 
A twenty-five hour course in programming for Whirlwind I was completed 
by the authors to partially fulfill this requirement. 

Once the capabilities of the computer were realized, it was 
necessary to decide what problem x-rould best prove the basic contentions. 
Probably the biggest single factor in this decision was time. A very 
complicated problem would require a great deal of time to program while 
a simpler problem would not test the authors' ability. 

A compromise was reached by deciding to program two problems. 

The first should be a relatively short and uncomplicated problem aimed at 
giving the authors experience with the computer and also making some 
contribution to naval architecture. The second problem should be of such 
size and complexity that its solution would more fully describe the 
abilities of the computer and the benefits of computer applications to 
naval architecture. 

The first problem programmed was the determination of the 
whirling frequencies of propeller shafts. The method of solution was 
developed by Professor F. M. Lewis. This problem was selected because 
the method of solution is an iteration process ideally suited to computer 
application. Furthermore, the problem was short and each step had been 
delineated by hand calculations performed by Professor Lewis. This last 
reason was very important because it enabled the authors to check their 
work and to ascertain the accuracy of the program. 
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The second problem programmed was the determination of stresses 



in a transverse web frame. A procedure similar to the one used in the 
Bureau of Ships was adapted to the computer. The numerous computations 
encountered in this problem made a computer solution very desirable. The 
calculations were made for a new destroyer type for which a partial check 
of the numerical results was available. Another important reason this 
problem was selected was to further the authors' knowledge in structural 
analysis of a ship. The course of instruction in naval architecture had 
not included the calculation of transverse strength. 



frequency of a propeller shaft system in which the excitation is induced 
by the propeller. The method of calculation involves a numerical solution 
of a fourth degree differential equation of the form 



PROPELLER SHAFT WHIRLING FREQUENCIES 



Description of the Problem 



Basically, the problem is to calculate the whirling or lateral 




h 2 



of which the general solution is of the form 



y = sin Xx + cos Xx + sinh Xx + cosh Xx 



where 




XL = p 
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let y* s. 9 and y" = Y. Applying boundary conditions to determine the 
constants of integration: at x = 0 



c -±L 

L 2 ~ 2 
* 2\ 



¥ 



Y 
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C. = 



C = 



2X 



7)(cos p-cosh p) -~sin p 



sinh p-sin p 

^(cos p-cosh p) +- 7 - sinh p 

2X_ 

sinh p-sin P 



at x = L 



„ V°s p cosh P-1, . 9 0 (co3 P 9 cosh ^ 

®L ’ si 4 p-slnp > * s:u * p-sta p 



Y _ Y(cos P sinh p-cosh p sin p) _ 2 Q\ sin p sinh P 
L ~ sinh p-sin p sinh p-sin p 



then writing 



«L ‘ - (S 0 S 1 * VV 

Y 1 ' - (S 1 Y 0 * V L V ‘W 
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From which 



S. 



1 



cos ft sinh p-sin P cosh [3 
sinh p-sin p 



q _ yj UUOH H~-L 

2 ~ p(sinh p-sin p) 



cos P cosh P-1 



2P sin p sinh P 




The above equations are applicable to one span of the shaft. 



In the solution of this problem these equations were used to determine the 
moment and angle of rotation at the end of one span and then carry these 
results over to the next span. The number of spans to be used is unspeci- 
fied, however, it is the opinion of the authors that four spans give 
sufficient accuracy. 



various values of p are assumed and the residual moment at the free end 
plotted versus the frequency. The value at which the end moment becomes 
zero is a resonant point or natural whirling frequency of the shafting 
system. 



Since 




In applying these equations it is useful to determine certain 



parameters in terms of the physical properties of the shaft system. These 
are given below and the notation is explained in Figure I. 
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Schematic Diagram of Shaft 
Figure I 



It can be shown from the equations 



previously developed that 



where 



and 



”bi = " l 
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m 



J = polar moment of inertia of propeller 
and entrained water. 

q = the order of vibration being investigated. 
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The number subscripts refer to spans as shown in Figure I, 
From the above it can be seen that 



”*L 2 “ " ^ s ll'"oi * °al S 31 ^ 

®12 ' ' (s uV + “oiV'W 



where 



c al ’ WH 



Here the second number subscripts refer to the span. It should 
be noted that 0^ is assumed as unity. 

For determining s^, s^, s ^ for different spans it is convenient 
to use the following: 

^2 = ^1 




where 



J dl 



E I 
1 1 

W L 3, 

1 1 



The procedure described is repeated for each end moment and 
angle of rotation until the residual end moment, m^ is determined. Then 
a new value of is assumed, and the process repeated until a series of 
end moments versus frequencies results are generated. By plotting these 
values, or by interpolation, the critical whirling frequency can be found. 
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General Discussion of Program 

The method developed for solution of the propeller shaft 
whirling frequencies is ideally suited to computer solution because it 
is essentially one of numerical iteration. Various values of are 
assumed and the residual moment calculated for any number of spans. Since 
the labor involved in programming is the same regardless of the number 
of times a process is repeated, good accuracy can be obtained with no extra 
effort by the programmer and small increase in cost of operation of the 
computer. 

The underlying guide in programming this problem was to make 
it as general as possible. This was done so that anyone wishing to use 
this method in the future would merely have to change a few numbers to 
adapt this solution to his particular problem. The quantities to be 
changed were grouped separately at the beginning of the program and 
involve such factors as the physical characteristics of the propeller 
and shafting, the number of spans to be considered and the order of the 
vibration being investigated. 

The first step in programming was to assign floating addresses 
or "flads" to registers which were to hold the numerical quantities and 
key instructions. These flads contain the numbers or instructions upon 
which mathematical operations are to be performed. The purpose of 
floating addresses is to simplify programming by enabling the programmer 
to refer to a flad rather than an absolute address each time an instruc- 
tion is written. The advantage of this system is clearly demonstrated 
in the programs contained in the appendix. The flads used in pro- 
gramming and the numbers they represent are given in Appendix B. 



In addition to storing the data and constants of any problem, 
registers must be assigned to store intermediate results and for temporary 
storage. 

The addresses of key instructions were assigned flads to aid in 
writing the program. For example, the flad, hi is the address of the first 
instruction in the program. After all the information has been fed into 
the computer the following instruction is used, "START AT hi" and the 
program begins 0 

After the flad assignments had been made, the next step was to 
draw a flow diagram such as shown in Figure III. The purpose of this 
diagram was to aid in the organization of the program and to provide an orderly 
procedure for the programmer to follow. It was a picture of the problem 
showing the sequence of operations to be performed. The programmer then 
wrote a short routine that would do_ the calculations of each section shown 
in the flow chart. The last step, then, was to combine these separate 
routines into' one program that would solve the entire problem. Of course 
considerable care must be exercised in writing and combining these routines 
since the results of one section will influence the calculations in the 
following sections. 

When the values of frequency and residual end moment had been 
generated in the arithmetic element, they were transferred to a magnetic 
tape via an output procedure. When the program had been completed these 
results were then transferred to a flexowriter and typed. The results were 
presented in two columns of frequency vs. end moment. Thus, it was possible 
to determine the frequency at which the moment becomes zero by either 
linear interpolation or a plot of the data. 



Figure HL FLOW CHART 



Pur /n Physical C oa/sta*ts aa/d /$, 



I W \ aTY~ 



store. //V f/ 



I Ca l 



STORL //V FZ 



\Cd -- WwlL 



store //v F 3 

31 



STORE AAV Ft 



Cj- = ^ (tj- 0 ) * 5 Toes /V Ff 



C F - ~'oeo)* Ybo', store aaa pi 
* ~ " 



c * =T 75 r 



STORE /A/ F 7 

X 



A* 


s ro^£ 


/A/ $/ 




^L_ 


sra^e 


JN & ^ 



ii 

> 

*> 


OUTPUT 


* 


<Pw=A*C w 


STO PE 


>N 3 3 


l 


? 

n 

& 

St 


STOKE 


AN g? 






5/ , S A) S 3 , STORE 
— . ilr ' 

Mu s - ( 5 / n*, - c A , s, e 0/ ) ’* 



a 






{ 5 £ dj>±. r 5 ^ | 

i _C^y J 


1 


sro/*£ //v ^ 




A = A A 


STORE AAA e 7 




/V?v //V & 6 ' 




OUTPUT 


4 ? 


TNCREA SE 


A 


sy Tiy o 


'if 


SELECT HE W 


WLUE OF $ 



17 



STRESS CALCULATIONS FOR A NON-SYMMETRICALLY LOADED XtfEB FRAME 
Description of the Problem 

A transverse web frame is a part of a system of frames and 
bulkheads that maintains the shape of the hull, A sketch of a web 
frame is given in Figure IVa, To do this it must provide support for 
longitudinal strength members and carry vertical shear loads from the 
decks and hull bottom into the sides which are the main vertical shear 
carrying members of the structure. The web frame is a continuous beam 
ring that encircles the hull girth in a transverse plane. All loads 
applied to the frame are assumed to lie within this plane, but may have 
any orientation within that plane. Under such loading a web frame must 
resist axial loads, shear loads and bending loads. Since web frames are 
generally very long with respect to depth, bending moments are by far the 
most important aspect to be considered in the solution. 

Due to the complex interaction between the webs, longitudinals, 
and bulkheads, the exact solution of the problem of web frames and the 
system of members that maintain the hull shape has defied rational engin- 
eering theory up to this time. Time limitations and the intent of this 
thesis prohibited an investigation of possible improvements in the method 
of the transverse strength calculations. It was decided to adapt to a 
computer solution a proven method. Although transverse strength calcu- 
lations are not always performed for modern ships, the analysis normally 
utilized, if performed, is a long, laborious method developed by J. Bruhn' J 
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Figure XVb Web Frame Showing Redundant Forces 
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Bruhn's method is an adaptation of Castigliano ' s theorem of 
least work which states that "The redundant reaction components of a 
statically indeterminant structure are such as to make the internal 
work a minimum. " ^ Determination of each of the redundant components 
requires that the differential of the total internal energy with respect 
to that redundant be equal to zero which is a minimum energy condition. 
Bruhn neglected the axial and shear components of energy and considered 
only the larger bending effects. Bruhn's basic energy equations for a 
Ttfeb frame with txro stanchions, as shown in Figure IVb, then become 




where m = m^ + p^ y + q^ x - m - m^ + R(x-x & ) + S(x-x^) 
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Bruhn's working equations become: 
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Gastigliano 1 s first theorem states that, "the deflection in a 
certain direction at a certain point of a statically determinate structure 
is equal to the partial derivative of the total external work or the total 
internal energy with respect to a load applied at this point in the 
direction of the deflection* " Conversion of a statically indeterminate 

structure into a statically determinate structure by means of cutting 
the structure and applying redundant forces as loads allows the application 
of Castigliano 1 s first theorem. The above equations then represent the 
deflection between the limits of the integration. Since the web frame is 
continuous, the deflections and rotations at the point of action of the 
redundants m^, p^ and q^ are equal to zero. For the redundant column 
forces R and S, the columns are considered incompressible and the vertical 
deflection between the column ends is equal to zero. This idea of 
deflections and rotations is used in the description and analysis of the 
problem in Appendix B. 

To determine the redundants, five simultaneous equations must 
be formed and solved. In a problem that contains more or fewer redundant 
supports such as stanchions or decks, the number of redundant equations 
required to solve the web frame would of necessity be changed to form 
one equation for each redundant force. 

After the equations had been developed and formed into a matrix, 
Crout’s method of matrix reduction was used to solve for the value of the 
redundants. These redundants were then combined with the applied loads 
and the moments recomputed. Utilizing the computed moments, the stresses 
due to bending were calculated. The direct forces were converted into axial 



22 



and shear stresses. These stresses and the bending stresses were then 
added algebraically to give the total stresses in the flange of the web 
frame and shell of the ship. 

General Discussion of Program 

The flow diagram contained in Figure V describes the sequence 
of calculations used in calculating the frame stresses. A detailed 
description of the instructions utilised in each sub-routine and the 
sequence control program are contained in Appendix B. 

Due to the normal manner of construction of ships, the web frames 
are constructed as a symmetrical ring frame of variable cross-sections. 

The loads, however, are very complex and for a general solution must be 
considered as non- symmetrical. Note that for this problem only the simple 
static loads are considered. 

The contours of the ship’s hull must be represented within the 
computer and, in order to find the neutral axis and to form the components 
of the loads, the slope of the contour must be computed at each station. 

The most obvious solution for the novice programmer is to form the equation 
of the contour and differentiate that equation to determine the slope. As 
novices this solution was attempted, but was soon abandoned due to the 
complexity of the program necessary to form the equation and the inherent 
error involved in first using approximate methods to form the equation 
and more approximate methods required to utilize the equation. The method 
finally selected was to choose a series of stations around the girth of 
the web and assume the web is represented by a series of short straight 
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beam segments. Stations were selected in such a manner that the angle 
change between segments was both small and nearly equal which, due to 
the variable curvature, makes each segment unequal in length. A station 
was necessary at the point of action of each redundant. Such a proced- 
ure follows a common philosophy of digital computers in that where 
complexity must be weighed against a greater number of calculations 
and simplicity, the choice must be made in favor of simplicity. The 
number cf stations required to properly approximate the DD hull was 
k9 which should be sufficient to approximate the frame of any vessel. 

Sin 9 and cos 9 were determined directly from the coordinates 
by using the tythagorean Theorem and basic trigonometric relationships. 

As = 

This was done for two reasons, first, only the sin and cos are used in 
the solution and secondly, the tan 9 becomes infinity when Ax = 0 

which is the case on ships with wall sides or tumble home. The compviter 
will not handle division by zero and the operations will be stopped ■when 
such a cas^ occurs. To approximate the functions of the tangent angle at 
each station, the functions of the chord angles were averaged. While 
sin 9^ + sin 9^/2 does not equal sin(9^ + 9^)/2, the error will be small 
as long as the change in 9 is small which is another example of using a 
simple idea many times to obtain accuracy and generality. 

Using the total water head -which in the present problem was 
the height of the main deck and the average of the y coordinates of the 






dx + dy sin 9 = 



^7 

As 



cos 9 = 



Ax 

As* 
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ends of the segments, the water pressure loads were determined and vertical 
and horizontal load components formed. 



f 

w 

V 



w 

H 

w 



(rgy. 2 ) 

- H- —5- 

= p ds 
w 

= f sin 0 
w 

= f cos 0. 
w 



Using the inputs of the web dimensions, the distance of the 
neutral axis (d^, ^ ) from the hull shell was calculated and the moment of 
inertia about the neutral axis determined ( 1^ ^ ) 

m t 

d N.A. = A t 

■Si.A. = I shell “ Wa. 

Knowing d M . and the (sin $) and (cos 6 ) the offsets of the 

IM • A m 

neutral axis from the shell could be determined* 



X N.A. * x s' d N.A. (sin ^ 
y M.A. - V d H.A. (cos 



Using the same procedure developed for the shell, the segment lengths 
and functions of the angle were formed. 
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One difficult engineering decision was the consideration of the 
sharp comer at the junction of the side shell and the deck edge. In 
normal practice the comer is bracketed and represents a point of complex 
discontinuity which nearly defies analytical analysis. After a number 
of methods were considered, it was decided to make the comer a "hard 
spot" by including the full bracket depth. The station spacing near the 
comer was reduced to minimize the effect of the comer. In all cases 
the tangent at the corner was taken as -1 which is equivalent to a slope 
angle of 135 degrees on the starboard side. This was done to form a 
smooth transition of the neutral axis around the comer. 

The moment due to applied loads and water loads was determined 

by m = H dy + Vdx. The moments due to the redundant forces were: 
n 

p^y + + R(x-x a ) + SCx-o^). 

The integrals of the moments were formed by the trapezoidal 
rule. While the trapezoidal rule implies linear segments between moment 
ordinates, this method was necessary to accomodate the random station 
spacing. The error resulting from the trapezoidal rule will be small as 
long as the station spacing is close and the moment curve is smooth. Both 
conditions are present in this problem. 

After solving for the redundants and recomputing the moments and 
axial loads to include all forces, the direct stresses and bending stresses 
were combined. The program as written rives the alternative of performing 
a check calculation on the set of frame scantlirrs or of altering the frame 
scantlings to approach a more efficient design. The stresses were compared 
with allowable stresses and the section altered by changing the web flange 
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Figure VI Arrangement of Matrix Formed by the Integration Sub-routine 
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■width by predetermined incremental amounts. If (^Tpiov” ^ > then the 
flange was reduced. If ( CT^ - Op < 0 then the flange was increased. 

The outputs of each cycle change may be any commuted value the 
programmer desires, but in the present problem the outputs were flange 
width, stress at the flange and at the shell, direct stress and shear 
stress at each point. Other methods of scantling change could have been 
utilized or the scantlings could have been altered by a complete change 
of scantlings around the girth. In the latter case each cycle would have 
required a decision by the programmer, a procedure that has considerable 
merit where design experience in the mind of the programmer can be better 
utilized. In the present problem it was decided to use the conditional 
response as another example of the digital computer’s capabilities. 

In the web frame problem the quantity of generated numbers 
necessary to produce a solution was so large that the auxiliary drum storage 
had to be utilized. Blocks of numbers were moved from the core memory onto 
the auxiliary drum to make room for calculations. Further, the size of 
the program was so large that optimum data output procedures could not be 
utilized without unwarranted complication of the program. The output data 
could have been typed or presented on films in columns of easily discernable 
numbers. In the present case the size of the program required output via 
an ''fp 1 ’ tape which gives five columns of numbers running consecutively 
across the lines and the address of the registers containing the numbers 
listed in the first column. Judiciously placed zeros provide separation 
between numbers of different meaning. 
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To make the program as general as possible sub-routines were 
utilized to perform the calculations. A full description of the use of a 
sub-routine is contained in Appendix A. The sub-routines used will perform 
the calculation of stresses for a web frame of any dimension and shape 
desired. However, the number of redundants of the frame cannot exceed 
five and the water pressure head must equal the height of the deck edge. 
Due to limitations of storage space, the s equence instructions used will 
permit the use of only k9 stations around the girth. All of the above 
limitations can be exceeded by minor changes in the program should the 
necessity arise. 
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RESULTS 



The numerical results obtained from the computer solution of 
problem one are contained in Table I. Numerical results were not obtained 
for problem two. 

More important than the numerical resui ts are the programs 
themselves. Problem one has been run and is known to be correct. The 
sub-routines of problem two, with the exception of the sub-routine for 
forming the stresses, were rim and found to check. When all of the sub- 
routines were run in conjunction with the sequence control program, several 
minor programming errors were found. 

The first program will determine the critical whirling frequency 
of any propeller shaft where no more than four spans need be considered. 

It can be adapted to another shaft merely by changing the physical con- 
stants of the shaft involved. 

The second program will, if the above mentioned errors are 
removed, determine bending, axial and shear stresses in a transverse 
web frame and, if desired, will adjust the size of the web frame until 
the computed stresses approach the allowable stresses. The program is 
general within the limitations stated below; the required inputs being 
the dimensions of the web frame, the location of the stanchions and 
coordinates or points on the transverse section, the weight and shear 
loads at these points and the hydrostatic design head above the keel. 

The following are the limitations of the program: 
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(1) the frame scantlings are symmetrical about the ship centerline 

(2) the points defining the transverse section must be symmetrical 
about the ship centerline and only 25 may be taken on one 
side ; 

(3) one of these points must be on the centerline keel, one 
on the deck centerline and one at each end of all the 
stanchions; 

(li) the hydrostatic pressure head must be equal to or greater 
than the height of the deck edge and hydrostatic pressure 
loads are not generated for the deck; 

(5) ship must be in upright position; 

(6) only one centerline or two off-center vertical stanchions 
may be used. 



Computer time for calculation of problem one was 10 minutes, 
but the typing time of the results was 20 minutes. Computer time for 
the solution of problem two was 10 minutes. 



Frequencies from 
hand calculations 



Frequencies from 
computer solutions 



+6 ORDER 



3.69 cps 


3.870 cps 


5.ia 


5.656 


8.20 


8.527 


11 .I 19 


11.762 


15 . 00 


15. 239 


21.00 


20.915 



-6 ORDER 



3.870 

5.656 

8.527 

Not Available 11.762 

15.239 

20.915 

+1 ORDER 

3.870 

5.656 

8.527 

Not Available 11.762 

15.239 

20.915 



Critical Whirling Frequencies for Propeller Shaft 

Table I 
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During the period that the thesis was being printed, the pro- 
gramming errors in the second program were remedied and the results listed 
in Table II were obtained,, Table II contains the values of the redundants 
and the stresses at the keel and the deck at centerline,, Due to the 
difference in station spacing between hand and machine calculations, the 
centerline points are the only coincident stations for the two methods 
of calculation. 







Hand Method 


Computer Method 


m 0 




+72,lhk -ft. 


+b6, 395 -ft. 


p o 




+20, 8ia 


+19, 805 


^0 




(assumed zero) 


+282 


R 




+ 1925 


+662 


S 




(assumed equal to R) 


+662 


11 

*!> 




-3880 # /sq.in. 


-5193 */ s C*i n » 


^flange 


>■ keel 


+30,750 Vsq.in. 


+ 18,970 Vsq.in 


^shell ' 


5 


-l5,UlO # /sq.in 


-16,586 # /sq.in 


' \ 




+i;70 ^/sq.in. 


-1,000 */sq.in. 


( deck 

u f [ centerline 


- 10 ,U 80 **/sq.in. 


+25,000 ^sq.in 


(T ^ 
s 




+3, 350 # /sq.in. 


-12,868 ^/sq.in 



Results of Web Frame Stress Calculations 



Table II 
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As shown in Table II there is little correlation between 
hand and machine calculated values. However, this lack of comparison is 
caused by two conditions. The most important condition is that the 
hand calculations have been found to be of doubtful validity. Several hand 
calculation errors were found and the scantlings used are not defined. 
Recent information indicates that the scantling values used in the 
Bureau of Ships hand calculations do not correspond with the scantling 
values presented in the contract plans used as a basis for machine 
calculations. 

The other condition that causes a small error is the fact 
that summation of vertical forces does not equal zero. The difference 
between weight forces and buoyant forces must be balanced by vertical 
shear forces applied along the vertical sides of the hull. This shear 
force was taken from the hand calculations and was included in the 
input values for the program. The machine computed buoyant forces were 
different from hand computed values due largely to the difference in 
station spacing and the method of averaging the water pressure. Therefore, 
the input values of shear did not balance vertical forces computed by 
the machine. It is now obvious that the shear force values must be 
computed by a sub-routine in the computer program. 

This lack of the means of proper comparison points out one 
other general rule of programming. The programmer should perform hand 
calculations for the problem being programmed using the method adapted 
to the computer. The authors decided that the use of existing hand 
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computations would be a time saving device ; however, as the results 
indicate, this was a poor decision and it is hoped that future programmers 
will benefit by the author's experience. 

With the exception of the shear error discussed above, the 
program does work and produces good results. With slight correction the 
web frame program is suited for research and engineering calculations. 

The time required for machine calculations was 25 seconds; 
however, program conversion required five minutes and typing time for 
results of calculations required twenty minutes. 

Time limitations did not permit the running of the comparison 
and alteration program for the web frame. 
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DISCUSSION OF RESULTS 

The subject of computer application to naval architectural 
problems is as general as the field of naval architecture itself. Each 
proposed application will require a separate analysis of feasibility. 
However, certain general conclusions and impressions can be gained from 
the two problem examples of this thesis. To each general conclusion there 
exists certain possible exceptions; however, these exceptions only tend 
to prove the contention that digital computers are not universal panaceas, 
but only local aids. 

The programming errors found in problem two are normal errors 
that occur in every program and do not reflect upon the method of approach. 
The nature of these errors are incorrect addresses, typographical mistakes 
and human errors. The correction of these difficulties is only a matter 
of time and patience. 

The results of the two programs indicate that digital computers 
are applicable to long and complicated naval architecture problems and 
that the numerical results are more than comparable with those of hand 
calculations. Of greater engineering importance than accuracy, was the 
great reduction jn calculation time obtained by computer usage. For the 
first problem the computer required 10 minutes as opposed to an estimated 
hand calculation time of thirty (30) hours and in the second problem a 
more defined superiority is estimated in that the commuter solution will 
require approximately ten (10) minutes versus an estimated sixty ( 60 ) 
hours by hand. However, this comparison does not depict the true time 
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picture and the long tedious hours of program writing and ever present 
"debugging" or program correction must be considered. Although the pro- 
gramming time varies greatly with the programmer's familiarity with the 
problem, computer and also, with the nature of the problem itself; in 
general the programming time and the time required for one hand calculation 
are roughly commensurate. The authors estimate that a practicing naval 
architect very familiar with the proposed problem and experienced in 
programming could program problem one in twenty (20) hours and problem 
two in forty (bo) hours. 

The "debugging" time for each problem would require another 
ten (lO) hours. The fact that the authors, unfamiliar with both the 
problem and computer, required six times that amount of time indicates that 
experienced personnel are a very important application factor. In the 
present case problem understanding and organization were responsible for 
the greatest loss of time. As a general statement it can be said that for 
an engineer to properly program a particular problem, his knowledge of 
that problem must be more thorough than that required for hand calculations. 
Not only must he know the individual steps in the problem, but also must 
predict the effect of each step upon the following step and the general 
value of the numbers generated. 

Also, the authors found that one programmer working alone is 
just as effective as several working together on the same problem. The 
reason for this is the necessity for a program to follow a single trend 
of thought and the exact accounting required. 
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To provide a more complete time description of nrogramning, 
the nature of the proposed problem should be analyzed. Problem solutions 
that involve repetition of the same equations i-rill tend to have a much 
shorter programming time. For instance the time required to program 
problem one for bO spans of shafting would be only slightly longer than 
that required to program 1| spans. Also, the program required to design 
the web frame by cycling and altering the scantlings is little longer 
than the program for check calculations. 

After the solution of one particular problem by a computer a 
general program will find further use in solving other similar problems. 

The programs written for this thesis are reasonably general and the 
limitations stated in the results are not of serious damage to the program. 
The limitation of problem one to four shaft spans in a reasonable engin- 
eering division in that propeller excited whirling vibration amplitudes 
beyond span four are so small that their effect on the end fixity of 
span foui is insignificant. In problem two the limitations were imposed 
to save storage space and reduce the complexity of the program. "This 
coincides with the authors' opinions that for engineering applications 
several simple programs are easier and faster to handle and program than 
one complex all purpose program. For instance, web frames are normally 
constructed symmetrically, 23 points should properly define most normal 
web shapes and to have a station at each major load or redundant is 
reasonable. The design hydrostatic head is normally above the ms, in deck 
and waterloads are normally considered a constant value along the deck. 
While the listed position of the ship is an important consideration. 
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modern naval architects consider the upright position for practical 
calculations and should the listed condition become necessary for re- 
search calculations, rotation formulas and non- symmetrical water loadings 
must be incorporated into the program. To change stanchions and add 
decks as redundant supports major alteration of integration procedures 
would be required. 

A general program is useful in design and research to compare 
and generate great quantities of data and to test the results of slight 
alterations upon a major problem. For instance, program one could analyze 
and compare the whirling frequencies of hundreds of ships in a relatively 
short time or the effects of propeller weight, fixity and blade number 
may be determined for numerous combinations of values. Problem two 
would find use in generating design schantlings for various load and size 
variations such as those found in classification society publications. 

As a sub-routine in conjunction with a more complex program the stiffness 
of web frames relative to transverse bulkheads could be determined by 
application of the theory of beams on elastic foundations. By reiteration 
of calculations under varying loads the variant or impactive loadings 
may be analyzed. Since the longest time consumed in a program is in the 
program writing, problem repetition increases the efficiency of the com- 
puter application. 

To be feasible for engineering calculations the computer must 
save time, money and increase the effectiveness of the engineers it aids. 
Therefore, the computer must accept numbers of engineering magnitude and 
quantity in a convenient form. The coded instructions must be straightforward 
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and have names descriptive of their use. Also, the generated results must 
be presented in such a manner that wide spread publication is possible. 
Whirlwind I and some commercially available machines will accept decimal 
digits of sufficient length and quantity to perform most practical 
engineering and research problems. Of extreme importance is the necessity 
of maintaining a large library of sub-routines capable of performing such 
calculations as sin ^9, etc. Frequently it is desirable to purchase 
the services of an experienced programmer to act as a consultant when 
first using a computer. 

In computer application there exists two basic schools of 
thought. One group thinks that problems should be programmed by experts 
who perform no other function than programming. In this case the engineer 
presents the problem to the programmer and at some later time is riven a 
solution. The other group states that a computer is best utilized by 
having engineers familiar with the problem effect the programming. While 
there exists a middle road somewhere between the two extremes, the awth-v 
are of the opinion that engineers should program the computer. This 
performs two useful functions for the profession as a whole. One, the 
programmer can make the necessary engineering decisions to adapt the 
problem more successfully and secondly, the engineer becomes more aware 
of possible computer applications in his field. 

The qualities required of a programmer are no more exacting 
than those required of a good naval architect. He must have almost 
infinite patience, imagination, ingenuity, good ability to organize 
problems and a good background in both cheoretical and numerical methods 
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of mathematics. The work is tedious and all instructions must be perfect 
since the machine will generate wrong answers as fast as correct ones. 
Digital computers are often referred to as "brains, " but nothing could be 
more incorrect in that thought is impossible. Actually a computer more 
nearly resembles a desk calculator or an abacus which operates at a hieh 
rate of speed. Imagination and ingenuity are so necessary to produce a 
workable, general program that many commercial computing companies obtain 
their programmers from chess clubs, cross-word puzzle enthusiasts and 
other hobbies where deep concentration and original thought are required. 

The future of digital computers in naval architecture appears 
very promising from both the research and practical viewpoints. In 
research they can be used for the rapid assimilation and processing of 
large quantities of technical data and for the solution of problems that 
require an inordinate amount of time to calculate by hand. Typical examples 
of the d ata processing type of problem would be a comparison of stability 
among several different vessels and other similar problems involving 
statistical analysis to establish design criteria. Since the same pro- 
gram may be used to solve a number of similar problems, design data in 
the form of charts and tables can be generated for the grillage beam 
analysis of decks and other complex problems. Such a procedure permits 
the determination and testing of important parameters for later applica- 
tion to specific designs. 

For specific design applications commuters may be used to effect 
large savings in time and labor. This is possible because the computer 
can approach more closely an optimum design with little increase in 
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computing time and expenditure of effort. In this category of appli- 
cation are midship section design, transverse framing design and flight 
deck bents. In addition, refinements intheory and in the extent of calcu- 
lations can also be utilized to a greater degree than practical by hand. 
Examples of these refinements are the use of grillage beam analysis for 
decks and bulkheads, the performance of stress calculations for a ship 
in drydock, design of propellers by circulation theory, and the investi- 
gation of submarine dynamical stability. In addition to design calcu- 
lations a digital computer can be applied very readily to scheduling 
and accounting problems. In shipbuilding this system can be used in 
determining the weight and center of gravity of a ship as it is being 
built. Another possibility is merchant ship ca.rao arrangement. 

Theses present an opportunity for wide and varied use of 
digital computers. In this application the computer should be used as 
an aid and not as the general subject. A student attempting a problem 
in which the magnitude of the calculations is the dominant factor should 
consider computer aid. A desireable arrangement is to have a team with 
one member learning the problem and another member learning the computer. 
In this manner both can derive maximum benefit from the thesis and a 
larger contribution to the profession can be made. Since computers are 
becoming more widely used, familiarity with their capabilities may be 



of considerable future value to the student 



CONCLUSIONS 



1. Digital computers are applicable to naval architecture problems. 

2. The average naval architect has sufficient background to program 
a digital computer. 

3. Some general requirements of a programmer are: 

a. They should have a complete engineering understanding 
of the problem to be programmed. 

b. Familiarity with numerical methods of mathematical 
processes is necessary. 

c. They should be amenable to tedious work. 

li» The problem to be programmed must be suitable to comcuters since 
it may take longer to program than to determine one solution by 
hand. 

5 . One programmer can work faster than two or more. 

6. In general the limitations imposed on a problem a re not serious and. 
are no more restrictive than those in hand calculations. 

7. Complete generality of a program is not always feasible. 

8. For engineering applications the following features of computers 
are necessary: 

a. Computer word-length should be able to handle engineering 
numbers. 

b. Library of sub-routines should be available. 

c. Computer should be easy to program. 
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RECOM-iENDATI QMS 



1. Digital computers should be used to solve naval architectural 
problems . 

2. Digital computers should be used in theses and research problems. 

3. Any future thesis should consider the computer as an aid in the 
solution of a specific problem and not applicability of computers 
to a general field. 

Iu For engineering applications the programmer should have; 

a. good engineering understanding of his field; 

b. knowledge of numerical methods; 

c. organizational qualities and a great deal of ingenuity; 

d. a working knoxfLedge of the computer utilized. 



Appendix A 

Supplementary Introduction 
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Supplementary Introduction 

A computer program is a series of logical instructions 
telling the computer what mathematical operations are to be performed 
on the data. Thus there are two groups of inputs, the data and the 
ins time tions. The data fed into the computer is stored in the core 
memory. Each number or piece of data is stored in a register or w cell ,r 
the location of which is designated by a number called the address of 
that register. After all the data is stored the instructions are 
stored in the unused portion of the core memory. The instructions are 
stored in sequence and will be performed in the order in which they 
are stored. No unused registers are permitted to interrupt the storage 
as these will stop the computer. The programmer may wish to know the 
address of key instructions such as the one that starts the computation 
process. He may assign an address to these instructions or allow the 
computer to do so, as explained later. 

The instructions given to the computer consist of two parts, 
the operation and address sections. The operation section is first and 
uses letters to tell the computer what mathematical operation is to be 
performed. The address section gives the address of the register contain- 
ing the data on which the operation is to be performed. The instructions 
recognized by Whirlwind I are listed in the glossary. 

When writing a program for a digital computer, there are several 
useful devices that simplify the procedure. The programming devices described 
below apply to Whirlwind I and may or may not apply to other computers. 

Probably the most useful device is a floating address. A 
floating address or flad is used to eliminate the necessity of referring 
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to the absolute register address in the core memory when writing the 
program. At the beginning of a program the inputs or key instructions 
on which mathematical operations are to be performed are given floating 
addresses. These addresses consist of a letter and a number, for 
example al, a2 etc. Any letter except 0 and L may be used and the sum 
of all the numbers behind the letters may not exceed 25>5>* In writing 
the program, then, the flad is used in the address section of the 
instruction to represent the absolute address. The computer will assign 
absolute addresses or register numbers, to these flads unless this is 
done by the programmer. For example, if al is the first flad in the prog- 
ram, the computer will assign it to register 32. If, however, the prog- 
rammer wants al in the register 130, he writes 130/al and the computer 
will place al in register 130. A table is usually produced by the 
computer showing absolute addresses assigned. It should be remembered 
that a flad can have only one value in any one program. 

Another useful device is preset parameters. These are numbers 
other than data that are used in the program. Preset parameters are 
frequently used to designate the number of times a computation is to be 
repeated. They differ from flads in that they may have any absolute 
value desired whereas flads must have the absolute address of a register. 
Preset parameters must be assigned a value before they are used and they 
will retain this value until it is specifically changed. Preset parameters 
are designated by two letters and a number. The first letter must be 
p, u or z and the second letter and the number raay be anything, except the 



letters o and 1 cannot be used 



In writing a program for a digital computer subroutines are 
very helpful. Basically there are Wo distinct types of subroutines which, 
for lack of a better name, shall be called internal and external. Internal 
subroutines are those performed automatically by the computer upon receipt 
of a single instruction. An example of an internal subroutine is the 
group of output instructions listed in the glossary. Each of these 
instructions will cause the computer to go through a certain predeter- 
mined procedure without any further information from the programmer. 

External subroutines are further divided into Wo categories; one category 
is designated library subroutines and are perfectly general for any prog- 
ram. The other external subroutines are groups of instructions that 
perform a repetitive section of a particular program. Library subroutines 
are groups of instructions usually written by an expert programmer for 
the purpose of performing a particular mathematical operation. They 
may be an integral part of the main program or they may merely serve as 
an auxiliary program. Examples of these external library subroutines 
are groups of instructions used to determine the square root of a number, 
trigonometric and hyperbolic functions of angles and the solutions to 
simultaneous equations. Most computer installations have libraries of 
subroutines available for programmers. These are written quite arbitrar- 
ily by using relative addresses and can be used with any main program. 

There are some rules that must be followed when using subroutines. 
The main program must be interrupted to call the subroutine into use. 

After the subroutine has performed its function, it must direct the con- 
trol element back to the proper instruction in the main program. An example 



of a subroutine follows: 



Major program 

ica a2 

f6, isp g3 — This instruction directs control to follow the 
instruction in the register designated by g3. 

SUB-ROUTINE 

g3, ita g3 plus 20 This instruction 

places the address contained in the 
control element in g3 plus 20 and adds 
one to that address. In this example 
the address contained in the control 
element is f6. 
iad t6 

more instructions 

tt ii 

ii it 

ii ii 

g3 plus 20, isp f6 plus 1 ■ — The address of 
the isp instruction was supplied by the 
ita instruction at the beginning of the 
subroutine . 

f6 plus 1, iad a3 

The authors do not wish to imply that a program should consist 
of a series of subroutines, but they are very useful in programs that 
perform similar operations many times throughout the problem. 



The numbers usually encountered in engineering are too large 
to be stored in one register which, it will be recalled, is l6-binaiy 
or about Ij decimal digits long. To accommodate these larger numbers 
a technique called floating point arithmatic is used. Two consecutive 
registers, giving 32 binary digits available for each number, are used 
to replace the single register. Some of these 32 digits are used to 
record the actual number and some to locate the decimal point. In this 
manner the range of significant digits can be increased by restricting 
the number of digits used to locate the decimal point. In the binary 
system the decimal point is written as some power of 2 and usually 6 
binary digits are used to express this power. Since the first digit 
of each register is used for the sign of the quantity in that register, 
in this case the sign of the second register is the sign of the exponent 
of 2. Twenty-four digits are left to express the significant number. This 
particular distribution of binary digits is referred to as the (2lu6) sys- 
tem and can handle numbers equivalent to eight decimal digits with expon- 
1 8 

ents up to 10 . Other systems are available and the size of the numbers 

in the problem determines which one should be used. 

In order to manipulate these two register numbers, an inter- 
pretive mode has been internally programmed to function as a part of 
Whirlwind I. All instructions that refer to double length registers 
are preceded by an i. If calculations or operations are to be performed 
on single length registers, the same instructions may be used without the 
i preceding. Hence there are two modes of calculation; interpretive 
for double register length numbers and whirlwind for single length register 



numbers. It is possible to pass from the whirlwind mode to the inter- 
pretive mode by use of the instruction IN and to reverse the direction 
by use of the instruction OUT. 



GLOSSARY OF INSTRUCTIONS 



TERM 



DEFINITION 



ica x 



ics x 



iad x 



clear MRA and add contents of 
register x. 

clear KRA and subtract contents 
of register x. 

add contents of register x to 
what is in MRA. 



isu x 



subtract contents of register x 
from what is in MRA. 



unr x 



idv x 



its x 



lex x 



isp x 



icp x 



multiply what is in MRA by the 
contents of register x and round 
off product to fifteen digits. 

divide contents of MRA by the 
contents of register x. 

transfer the contents of the KRA 
to register x losing what was 
previously in register x. The 
MRA remains unchanged. 

exchange the contents of the MRA 
with the contents of register x. 

take the next instruction from 
register x. Does not affect the 
contents of the MRA© An uncon- 
ditional instruction used to break 
the sequence of operations. 

if the contents of the MRA are 
negative, take the next instruc- 
tion from register x and continue 
from there. If the contents of 
the MRA are positive, take the 
next instruction in sequence. A 
conditional instruction used to 
break the sequence of operations. 
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Glossary of Instructions (Continued) 



TERM 


DEFINITION 


isc x 


select counter number x. Without 
this instruction counter zero 
will be used whenever +C appears. 


icr x 


cycle reset. Sets index register 
of counter to zero and criterion 
register to x. 


ict x 


cycle transfer. Increases con- 
tents of index register by one. 

If contents of index register 
greater or equal to contents of 
criterion register, set index 
register to zero and do next 
instruction in sequence. If con- 
tents of index register less than 
contents of criterion register, 
take next instruction from register x 


ici x 


increase the contents of index 
register by number x. 


icd x 


decrease the contents of index 
register by number x. 


iti x 


transfer the right eleven digits 
of the index register into the 
right eleven digits of register x. 


iat x 
ita x 


add contents of index register to 
the contents of register x and 
store the result in the index 
register and register x. 

replace the address section of the 
instruction in register x with the 
address that is one more than the 
address of the register containing 
the last isp or icp if the contents 
of the MRA are negative. 



ita x 



Glossaiy of Instructions (Continued) 



TERM 


DEFINITION 


iTOA 


record the contents of MRA on 
direct printer (typewriter). 


iMOA 


record the contents of MRA on 
magnetic tape for delayed printing 


iSOA 


record the contents of MRA on 
oscilloscope for photographing. 


iFOR 


this instruction provides an 
automatic device for obtaining a 
suitable layout of output data 
in columns, lines or blocks. 



Note: The above output instructions are usually followed by a series 

of letters and numbers that indicate the desired form and 
arrangement of the output. 



Appendix B 
Details of Procedure 



Details of Program for Determining the Whirling; 
Frequencies of Propeller Shafts 



52 



The first step in programming this problem is to put the 
numerical values of various quantities into the computer. This is done 
by the use of floating addresses, or "flads." The proper use and func- 
tion of flads is explained in Appendix A. 

The following table lists the flads, the quantity and the 
values used. Figure JL defines the symbols and nomenclature. 



al, 

a2. 



bl, 

b2. 



cl, 

c2. 



dl. 



h 



77" - length of tail section. 

513.6" - length of inboard sections. 



l 2 = 576.0" 
L 3 = 510u6" 



Vo 


= 17.U A 


10 10 


hh 


= 13.7 X 


10 10 


E I 
2 2 


= I2,h * 


10 10 


E I 
3 3 


— 12 ,h x 


10 10 


W 

P 


= U7,8i5 


lb - 




= .281; lb/in 3 



lb-in - E times moment of inertia of 
tail section* 

2 

lb- in - E times moment of inertia of 
inboard sections* 

lb-in 2 
lb- in 2 

weight of propeller 

- specific weight of material used in 
inboard shaft sections. 



^2 

*>3 

*1 



.281; Voj lx? 

,28b lb/in 3 
2 

193«5 in - average area of shaft cross section 
for each span* 



A 2 = 17ii.5 in 2 
A 3 = 17U.5 in 2 
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d2, 

d3, 

e3. 



el;, 

e6. 



3 - number of shaft spans used. 

135 - number of increments of p desired to 

be used. p is increased by increments 
of 2 degrees. 



<M> q - 6 



(i -i) 

k 2 q' 

( 1 - 1 ) 

V 2 q' 



q is the order of frequency being 
_ ^ investigated. Plus sign is used 



for backward precession and the minus 
sign for forward precession. 



q = 1 



J = 211.0 y. 10 - polar moment of inertia of 

propeller and entrained water. 

= 90 - initial starting value of ,8^, in degrees, 



The above listed flads are the only ones that must be changed 
in order to adapt this program to another problem. In order to accomplish 
this change it is merely necessary to add a tape with the correct values 
at the end of the program. 

The following table lists flads that contain constants which 
will remain unchanged for every problem: 



el. 


TT = 


e2. 


180 


e5. 


g = 


®7, 


Pl" 


e8. 


(— ) J 
'l 80 ' 


e9. 


TT 

(l8o)‘ 


elO, 


2 


ell. 


1 


el2. 


3 


el3. 


h 



g = 386 in-sec* 



(p^ + increment) 



A block of registers consisting of f and g flads has been 
set aside as intermediate storage. This is necessary to enable the 
computer to calculate quantities, store them temporarily and use them 
later in the program. Another group of registers consisting of t and u 
flads are set aside for temporary storage. These quantities are used 
once and then discarded. 

The last group of flads used is the h and m group. These are 
used to control the program by assigning key instructions to them. 

The outputs from this program will consist of a series of 
frequencies and corresponding end moments. These values will be arranged 
in two columns with the frequency corresponding to the moment that 
appears on the same line. It is very important to include the sign of 
the moment, for we are more interested in the point where the moment 
passes through zero than the actual values obtained. 

In performing certain routine arithmetical operations, such 
as finding the square root, use was made of the sub-routine library in 
the computer laboratory. These sub-routines are very reliable and save 
considerable time when programming. 

The actual orders used for programming this problem follow. 

The short paragraph at the beginning of each section explains briefly 



the function of the orders in that section 
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These three orders place the sub-routines for square root, 
sine-cosine and sinh-cosh at the beginning of the program. This is 
veiy desirable as it eliminates any confusion with the main program. 

Order 

51, LSR, FU8 

52, LSR, FUiia 

53, LSR, FU£ 

Put in flad assignments al through glO as previously noted. 

These orders are assigning values of 3 , 2, 3 in the flads 
Ml, M2, M3, respectively, the function of these flads is to set the 
cycle counter criterion registers to a predetermined value. The opera- 
tion is done in WW1 since the numbers are single register numbers. 



The last order is 


used to 


enter the CS 


mode. 




Order 




Input to MRA 


In MRA 


Transfer 
from MRA 


hi, ca. 


d2 


3 


3 




td, 


ml 




3 


3 


td. 


m3 




3 


3 


su. 


ell 


-1 


2 




td, 

IN 


m2 




2 


2 



This group of orders will calculate the weights of each shaft 
span and store the results in fl. The first order says that the opera- 
tion is to be repeated three times. The +C at the end of the next 



four orders is to permit the address section of the instruction to be 
increased by two registers on each cycle thus allowing three different 
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weights to be calculated by the same set of orders. The last instruc- 
tion is necessary to return control to the first instruction at the end 
of each calculation, to increase the address of each instruction by two 
registers after each operation, and to pass control to the next order 
when the operation has been performed the number of times specified by 
the first order. 





Order 


Input to MRA 


In MRA 


Transfer 












from MRA 


ml. 


icr 










h2. 


ica. 


c2+c 


A 


h 






ilTlTj 


a2+c 


* L 1 


V y i 






imrj 


dl+c 


X ^ 


V V \ 






its. 


fl+c 




vr 


If 




y 




1 


1 



This group of orders will calculate and store the results 
in f2. The cycle will be performed three times since the order ml, icr 
is still effective. 



Order 


Input to MRA 


In MRA 


Transfer 
from MRA 


h3, ica, b2+c 




Vi 

E„ I_ 




idv, a2+c 


+ h 


1 1 
h 

E I 




its, f2+c 




1 1 
h 


C 

a 



ict, h3 
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This group of orders will calculate for each span and 
store the results in f3. The cycle will be repeated three times since 
the order ml, icr is still effective. 



Order 



hi*, ica, f2+c 



idv, a2+c 



idv, a2+c 



idv, fl+c 



its, f3+c 



Input to MRA 

hi 

h 



In MRA 

hh 

h 

hi 

5 



E I 
11 

5 

hi 

hi 3 



Vl 

hi 3 



Transfer 
from MRA 



C 



d 



ict, hi* 



2 

This group of orders will calculate the constant n/(l8o) 
and store the result in e9. This operation will only be performed once 
as no cycle is indicated. 
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Order 


Input to MRA 


In MRA 


Transfer 
from MRA 


ica, el 


TT 


tt 




idv, e2 


— 180 


tr /180 




idv, e2 


-r- 180 


rr/(l80) 2 




its, e9 




tr/(l80) 2 


ir/(l80) 2 



This group of orders will calculate the constant (ff/l8o)^ and 
store the result in e8. The device used here is that when a number is 
transferred from the MRA to a register the MRA retains the number. The 
fourth power is formed by multiplying the number in the MRA by the number 
that was transferred which is the same as squaring the number in the 
MRA, in this case (tt/18o) 2 „ 



Order 


Input to MRA 


In MRA 


Transfer 
from MRA 


imr, el 


X TT 


(tt/i8o) 2 




its, e8 




(tt/i8o) 2 


(tr/180) 2 


imr, e8 


X (rr/180) 2 


(tt/180)^ 




its, e8 




(n/l8o)^ 


(tt/180)^ 



This group of orders will calculate C and store the result 

w 



in fij. 



Order 



Transfer 
from MRA 



Input to MRA 



In MRA. 



ica, f3 



imr, cl 



imr, al 



imr, al 



imr, al 



idv, bl 



imr, e8 



its, fli 



Vl 



a w 
p 



X L„ 



x L, 



X L 



‘0 



-Vo 



/ rr J 

* 'iso' 



5 A 

"A 



Wi 

Vi 



w pWi 

W A 3 

W L^E I 

P Q i i 

Vl 

Vfol 

Vi 

W L^E I_ 



W. 



P 0 1 1 

'aVo 



w 



^oh 1 ! , " 'll 

VlVo 180 

VfcA , « )h 

'l8o' 



m 



w 



This group of orders will calculate and store the result in 
f6. In this operation the square root sub-routine is used. The third 
order is used to enter this sub-routine, which is in the Whirlwind mode, 
requiring the IN order to return to the CS mode once the square root has 



been calculated 
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Order 



Input to MRA 



In MRA 



Transfer 
from MRA 



ica, e£ 
imr, f3 

sp, si 



IN 



imr, e9 



idv, elO 



its, f6 



g 



V i 
“A 3 



g 



X (l8o) 2 



-j- 2 



Vl 



Vl£ 

W Ip 
1 1 



¥ii 

Vi 



n 



|E I e 
1 IT 



(iSO)-/ WiI 3 



2(l8o)‘ 



TT 



hh* 



m u i l i 



2(180)2 /W 1 i3 



This group of orders will calculate and store the result 
in f7. The operation is performed twice, since m2 is 2, because this 
constant is used to determine and If more spans were used this 

operation would be performed a number of times equal to two less than 
the number of spans. The square root sub-routine is used twice in 
succession here to form the fourth root of a number. 
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Order 


Input to MRA 




In MRA 


Transfer 
from KRA 


m2, icr 










h 5, ica, f3 


¥1 

Vi 








idv, f3+2+c 


. Vs 

Ml 3 
2 2 


w i 




sp, si 










sp, si 










IN 




h 


w 2 L l 






\l 






its, f7+c 




k 


E L I 1 W 2 1 2 3 


o 




\l 


\ L l E 2 X 2 


c p 



ict, h$ 

This group of orders will calculate C^, and store the result 
in f5. This operation is to be performed three times as indicated by 
the first order* The flad ml was not used here because the number of 
times this operation is performed depends on the number of vibrational 
orders being investigated and is independent of the number of spans* 
Usually three orders of vibration will be investigated, (synchronous, 
forward and backward precession) . The flad h6 is the beginning of 
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the first loop of the program. This means that when the program has 
been min through once it will automatical l y come back to h6 and do the 
entire program again using the second value of e3. This loop will be 
performed three times as indicated by the first instruction. 



Order Input to MRA 



In MRA Transfer 

fron MRA 



icr, 3 



h6, ica, f3 



imr, eJj. 



imr, al 



idv, bl 



imr, e3+c 



imr, e8 



Vi 

Ta7 t3 



Vi 

W L 8 
11 

w 

T*r t.3 



X. L 



0 



E i “S JL n 

Vi 



tE-I. 
0 0 



W L o 

E oWi 






El , 

— ■ 1 ° 3 (i + J) 

E oWi q 



r) 



A Vi®! 



L 1_ JL^ _ , 

Jr . 1 — Q_ (i + 1.) (-2- 

wA 2 q 180 



( 1 t l )e J!_)U 

E 0 I q W 1 L^ 2 <1 180 



C. 

3 



its, f5> 
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This group of orders assigns a value to the flad ml*. The value 
used is that contained in d3 which for this problem was 135. This number 
was selected because we wanted p to vary from 90° to 360 ° in as small an 
increment as possible to increase the accuracy. Due to limitations on 
the magnetic tape output the smallest increment feasible was 2°. The 
first instruction is used to select the Whirlwind mode since the trans- 
fer of digits as indicated by the third instruction cannot be performed 
in the cs mode. The fourth order is used to return to the cs mode 
before proceeding with the remainder of the program. 

The last instruction is used to select a new cycle counter. 

This instruction permits a loop within a loop to be performed. In the 
preceeding group of instructions counter zero was used and is stepped 
once each time the sequence is performed, but a different counter is 
used to allow the following sequence to be performed, and counted, a 
specified number of times without interfering with the previous counter. 



Order 


Input to MRA 


In MRA 


Transfer 
from MRA 


OUT 


ca, d3 


135 


135 




td, ml* 
IN 




135 


135 



is, cl 

This group of orders will calculate the frequency and print 
the result on the magnetic tape output device. This operation will be 



6U 



performed 135 times, indicated by the first order, for each cycle of 

2 

the main program. The second through the fifth orders calculate p 

and p^ and store them in gl and g2 respectively. The frequency is then 

2 

found by multiplying C by P . The last instruction is an order to the 
computer selecting the output device and giving the information necessary 
to print the frequencies in a column, the number of digits desired, etc. 



Order 




Input to MRA 


In MRA 


Transfer 
from MRA 


mil, icr 










h7, ica. 


e6 


90° 


xo 

o 

o 




imr. 


e6 


* 90° 


(90°) 2 




its. 


gl 




(90°) 2 


<>2 


imr. 


gl 


* (90°) 2 


(90°)k 




its. 


g2 




(9 O 0 ) 1 * 


1 


ica. 


f6 


tt l E l 1 ! g 


« I s lh s 




2 (l 80 ) | ¥ 1 L^ 


2(180) 2 J 





imr, gl 

imoap 123lu56ssss 



X (90 °f 



n(9Q) 2 hh 8 

2(180) 2 / ¥ 1 L^ 



This group of orders will calculate the moment at the inboard 
end of the tail shaft, m^, and store the results in g5. This operation 
will also be performed 135 times for each cycle of the main program. 
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$ and ji J9 both functions of p are computed first and stored in g3 and gl| 
respectively. Using the formula in the flow-chart they are combined to 
calculate m^. ^t should be noted that this and the preceeding calcula- 
tions do not require more than two double length registers to store the 
results. In the frequency calculations the results are printed out as 
soon as they are obtained. Here the results stored in g5 are used immedi- 
ately in the subsequent calculations after which they are no longer desired. 
Therefore, the m^ calculated is lost as each succeeding moment is determined. 



Order 



Input to MRA 



In MRA 



Transfer 
from MRA 



ica fli 



imr g2 



its g3 



w lJel , 


¥ L 8 E T 


p 0 1 1 / n vJU 
W L 8 E I 180 


p 0 1 1 

W Pe I 


110 0 


110 0 




w lJei 

p 0 11 



180 



w iHVo 180 



w 



hhh /Ttp k 






ica 



imr e6 



W lWo 

f5 - 1)^)'* 

WiH q 

,1 i K w xii 0 

E I W i3 l 2 W p 

Vo"! 1 ! 



A 



w 



E I W l? 2 q l8 °" 
0 0 1 1 



gk 



,1 lw n xU R , 

Eq^lP 1 180 



its 



66 



Order 



isu git 



its t2 



Input to MRA 



In MRA 



(continued) 



ica 


g3 


r TiJ 


w 


idv 


el2 


-5-3 










3 


its 


tl 




/ 

r w 








3 


idv 


el3 


fit 


'w 

12 


imr 


gi* 


* 


12 


iad 


lh9 


+ 1 




isu 


tl 


t 

r w 

“ 3 








w 



1 + 



1 + 



J 



Transfer 
from MRA 



_1L 



12 



- 



4 



w 



12 



1 4 . 

5 J 1+ 12 



K j 
T" *1 



ica tl 



imr git 



i 



w 



3 

* h 
-h 



i 



w 



3 

'jA j - <t> ■ 



isu git 



Order 
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Input to MRA In MRA Transfer 

(continued) 



from MRA 



isu g3 



- A 



w 



w J_ - j&T - 



J 'w 



idv t2 ' -s- 1 + 



A. 



w 



12 3 






n VJ os. i 

1 + 12 3 " 



imr bl 



*Vo 



E_I_ 
0 0 



^w^J v v 
3 ~ h ~ K 



1 + 



'w 

12 3 " 






idv al -r L 



0 



its g5 



Vo 



- / T 

'w 



1 + 



W y 

12 “ 3 " *«U 



Vo 


1 

1 

I 


L o 


^w^ J ^w / 

+ 12 - 3 _ 6J 



m, 



01 



This group of orders will first convert the value of j3 currently 
stored in e6 from degrees into radians and store the result in e7« The 
last two orders will place the initial value of which is assumed as 

unity, in g9. 
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Order 


Input to MRA 


In MRA 


Transfer 
from MRA 


ica 


e6 


P 


P 




imr 


el 


TI 


"P 




idv 


e2 


180 


JIL 

180 




its 


el 




JIL 

180 


JIL 

l8o 


ica 


lh9 


1 


1 




its 


g9 




1 


O 

CD 



This group of orders will calculate various parameters, such 
as cos (3, sin [3 and store them temporarily foi future use in the calcu- 
lation of end moments. The first instruction selects a third counter 
to be used for this and subsequent calculations. This is done for the 
same reason as before, to permit this sequence to cycle three times, 
once for each inboard span, without affecting the cycle count of the 
previous loop. The second instruction sets the number of times the loop 
is to cycle. This number equals the number of spans since m3 takes its 
valud from d3 which is the flad assigned for the number of spans. The 
manner in which the subroutines are utilized is shown here and their 
advantage becomes evident. 
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Order 



Input to KRA 



In MRA Transfer 

from MRA 



isc 


2 








m3, icr 










h8, ica 


el 


Pi 


Pi 




isp 


s2 




COS 




its 


ul 




COS p^ 


cos P^ 


ica 


el 


Pi 


Pi 




isp 


ls2 




sin p^ 




its 


u2 




sin P 1 


sin P i 


ica 


el 


Pi 


^1 




isp 


s3 




sinh p^ 




its 


u3 




sinh p^ 


sinh p^ 


ica 


el 




Pi 




isp 


5>s3 




cosh p^ 




its 


u h 




cosh p^ 


cosh p^ 


ica 


u3 


sinh p^ 


sinh p 




isu 


u2 


-sin p^ 


sinh p^-sin 8^ 




its 


u3 




sinh p^-sin p^ 


sinh p^-sin p^ 


ica 


ul 


cos P^ 


cos p^ 




imr 


U3 


X. sinh 


cos p^ sinh p^ 




its 


u6 




cos P^ sinh p^ 


cos p^ sinh p^ 
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This group of orders will calculate the parameters si, s2, s3, 
and store the results in g6, g7, and g8 respectively. This operation will 
be repeated three times, once for each span, since the first two orders 
of the previous section are still in effect. 



Order 


Input to MRA 


In 15RA 


Transfer 










from MRA 


ica 


u2 


sin 


sin P x 




imr 


uh 


* cosh p^ 


sin p^ cosh p 




isu 


u6 


-cos p^ sinh p^ 


sin p^ cosh B^-cos p^ 


sinh p 








sin p_ cosh p -cos P„ 


sinh Pt 


idv 




4- sinh p^-sin 


1 11 
sinh 8^-sin p^ 


1 








sin p. cosh p -cos P 


sinh P„ 


its 


g6 




1 ] 11_ 

sinh p^-sin p^ 


L \ 


ica 


u3 


sinh P^ 


sinh p^ 




isu 


u2 


-sin P 1 


sinh p^-sin p^ 




imr 


el 


X P x 


P^(sinh p^ -sin p ) 




its 


u5 




P^(sinh p^-sin P 1 > 


P 1 (sinh P^ -sin p 1 ) 


ics 


ul 


-COS P-j^ 


-cos p^ 




imr 


uit 


X cosh p 


-cos p cosh p 
*1 *1 




iad 


lh9 


+1 


1-cos p^ cosh p^ 










1-cos p n cosh p_ 




idv 


u5 


-rp^(sinh p -sin p^ 


P-, (sinh p -sin p ) 





Order 



In MRA 



Transfer 
from MRA 



Input to MRA 

(continued) 



its 


g7 




ica 


u3 


sinh p^ 


isu 


u2 


-sin p 1 


its 


u5 




ica 


elO 


2 


i inr 


el 


,P 1 


imr 


u2 


X sin p 2 


imr 


u3 


X sinh p^ 


idv 


u5 


-rsinh p^-sin 



its g8 



1-cos p cosh p 

-i- -L c* 

P^(sinh P-j-sin p^) 2 

sinh p^ 

sinh p^-sin p^ 

sinh p^-sin p^ sinh p^-sin p 
2 

2 h 

sin ^ 

2P^ sin p^ sinh p^ 

2p^ sin p^ sinh p^ 
sinh p^-sin p^ 

2P n sin p^ sinh p 

J- X J- Q 

sinh p^-sin p^ 3 



This group of orders will calculate the end moments and the 
angle of rotation at the inboard end of the shaft selected. The previous 
results stored in g5, namely m^, and previously calculated values of 
9 qP si, s2, s3 are combined into the formulas shown in the flow chart 
to determine the values at the end of the second section, ^is 

process is repeated three times to produce the moment we are 

looking for. 0^ is also determined although it is not utilized. 
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The method of saving registers that is used in this section should be 
noted. is in g!?. As each successive moment is calculated it too 

is put in g5« The preceeding moment is lost by this process, but it 
is of no importance since only the final end moment is required. 

The last six orders have a different function. The first five 
of these are utilized to calculate the value of for the next span 

which is used in the calculation of the moment and si, s2, s3. The 
relation of f3 and (3^ is shown in the flow chart. 

The last order is used to cycle the program back to h8. ict 
causes the computer to return to h8 and at the same time increases the 
index of the cycle counter by one. When the index equals the number 
previously set in the criterion by an icr instruction ict will reset 
the index to +0 and pass control to the instruction immediately following. 



Order Input to MRA In MRA Transfer 

from MRA 



ica 


g5 


m 01 


m 01 




its 


gio 




ra oi 


m oi 


ics 


g5 


H 

e° 

1 


1 

o 3 

H 




imr 


g6 


* s i 


■"oi s i 




its 


ul 




- m 01 S l 


-”oi s i 


ics 


f2+c 


-c 

a 


-c 

a 




imr 


g8 


XS 3 


- C a S 3 




imr 


g? 


* e oi 


- C a S 3 S 01 
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Order 

iad ul 
its gf> 
ics glO 
imr g7 

idv f2+c 
isu g6 

its g 9 

ica f7+c 
imr e6 
imr el 

idv e2 

its e7 

ict h8 



Input to KRA 


In MRA 


Transfer 
from MRA 


(continued) 




•"oi s i 


“ C a S 3 6 0l“ m 01 S l 

“ C a S 3 G 0r m 01 S l 


*12 


" m oi 


~ m oi 




* S 2 


“ m 01 S 2 




-4-C 

a 


" m 01 S 2 

C 

a 




- s i 


- m 01 S 2 - 
c “ s i 

a 






■ m 01 S 2 , 
c ** S 1 

a 


9 12 


c p 


°p 




*p 


V 




>ur 


"V 




•I- 

H 

CO 

O 


355 V 






iso y 


^2 



This group of orders will print out the contents of g5, which 
are the moments at the inboard end of the shaft* The first instruction 



7k* 



is used to return the counting of cycles back to the 135. This is 
done since we want 135 end moments printed, or one for each B value 
used. 



Order Input to MRA 

isc 1 
ica g5 

imoa+p 123U.56 10 ^c 



In I--RA Transfer 

from MRA 



m 3k 



This group of orders does several things which will be dis- 
cussed in the order of their occurance. 

The first three orders are used to increase the value of P 
in e6 by 2° before the next cycle starts. This operation is performed 
135 times also, the p starts at 90° and goes to 360 ° in steps of 2°. 

The fourth order is used to break the normal sequence and 
add a correction before returning to the main program, a process known 
as patching. The correction added is given below. The first instruc- 
tion is merely a new flad used to get the number +1.0 in the computer. 
The second instruction increases the index of counter 1 by one and per- 
forms the same function as explained above, namely, permit this part of 
the program to cycle 135 times. The third and fourth instructions are 
merely to print a series of zeros in the print out to separate the 
results into blocks of 135 lines of frequency vs. moment for each order 
being investigated. The last instruction is to pass control to the 
next patch. 
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Order Input to MRA In MRA Transfer 

from MRA 



ica 


e6 


0 


iad 


elO 


+2 


its 


e6 




isp 


3h9 




isp 


hl2 




ict 


h6 





h9, i ST0P 

START AT hi 
lh 91 + 1.0 
ict h? 

ica lh9 +1.0 

imoa+123liCc 

isp-2h9 



0 

p+2 

P+2 p+2 



+ 1.0 



This group of orders will do several things. First it passes 
control to the flad hl2 which introduces the subroutine of resetting 
the value of p to 90° at the completion of the program for one order. 
Then counter zero is selected to begin calculations for the next order. 
The instruction OUT is used to enter the Whirlwind mode and the next 
three instructions are used to increase the address of l|h6 by two. This 
is normally done automatically by the use of counters and +<; at the end 
of each instruction. However, difficulty was encountered and this 
"’manual'’ stepping was used to save time. The last two instructions are 
used to return to the CS mode and the main program at one instruction 
before h9. 
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Order 


Input to MRA 


In MRA 


Transfer 
from MRA 


hll, +90,0 
•2h9/isp hl2 


hl2, ica hll 


+90,0 


+90,0 




its e6 




+90,0 


+ 90,0 


isc 0 
OUT 


caidi6 


28 $ 


285 




ade 11 


+1 


286 




ade 11 


+1 


287 




tsldi6 




287 


Uh6 



IN 

isp h9-l 



The main program is reentered as explained above and the last 
three instructions listed below are performed. The first increases the 
index of counter zero by one and when the entire program has been repeated 
three times passes control to the second instruction which stops the 
computer. 

The last instruction is used to start the calculations at the 
proper place. It is always the last instruction fed into the computer 
since the computer scans the entire program, locates the flads and stores 
the information before it starts. If this instruction were anywhere 
else the computer would start computing before all the information was 
fed into it and properly stored. 



Order 



Input to MRA 



In MRA 



Transfer 
from MRA 



ict h6 
h9, i STOP 
START AT hi 
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PROGRAM OF STRESS CALCULATIONS OF WON- SYM? TETRI CALLY 
LOADED TRANSVERSE WEB FRAME 

The working program of the web frame stress calculations will 
consist of a series of sub-routines joined together by a sequence control 
program (see Figure V). To facilitate programming and program alteration, 
the instructions and input data were placed in four separate sections 
and these sections were represented by the flexowriter tape number that 
contained the instructions on data. Scantling data and other information 
relating to a particular problem were placed in tape 359-217-100, The 
sub-routines were placed in 359-217-101, and, the sequence control pro- 
grams were placed in 359-217-221 and 359-217-222. Calculations for a 
different though similar ship will require the complete changing of 
tape 359-217-100. Tape 359-217-101 is general for similar problems. Tapes 
359-217-221 and 222 determine the nature of the calculations, that is 
whether a check calculation or a design calculation will be calculated. 

In the instructions for the input data in tape 359-217-100 abso- 
lute addresses were assigned to flads. The absolute addresses were necessary 
for block in, bl out procedures used in the auxiliary storage drum. The 
use of flads is still retained due to the ease of writing and categorizing 
letters. The flad and absolute address assignments for storage of 
generated data were also placed in tape 359-217-101. 

The following list of data are the specific input data for a 
DP type vessel. Due to the quantity of numbers in each ategory, only a 
sample of the total group is presented below. 
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fc Tape 359-217-100 

(2U,6) 


32 I xl, +00.00' 
+3,266' 


Assigns the first value of the x-coordinate of 
the shell to register 32 which is the first 


+6.508' 


register in the case memory. Note that the 


+9,917 


coordinates are in feet measured from the origin 
at the keel centerline. 



Total of 25 double length numbers which occupy 50 registers. 



+0.000 


Value at the main deck centerline. 


82 1 yl, +00.00' 


First value of y-coordinates of the shell. 



25 numbers 



+26.33 


y-coordinate of main deck centerline. 


332|bl, +2lj..00 


Value of web depth (d ) at the origin measured 

w 

in inches. First value stored in register 332 « 



25 numbers 



+12.00 

282 | b2, +0.3125 


Value at main-deck centerline 0 

Thickness of web (t ) measured in inches. 

w 


25 numbers 




+0.2500 
ii82 1 b3, +12.0 


Width of flange (w^) measured in inches. 


25 numbers 




1.0 

U32|bli, +0.500 


Flange thickness (t^) measured in inches 0 


25 numbers 




+0.250 
382 |b5, +0,7500 


Shell thickness (t ) in inches. 

s 



25 numbers 

+0,375 
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632 e5, -97k .0 



50 numbers 



-97it.O 

932 e2, +25.517 



978 el, +8*0 



pal = 25 



pa2 = 18 

pa3 = U 



Weight and shear loads at each station. The 
weights are determined in averaging the weight of 
vessel and shear for one foot of length. The length 
per foot times the number of feet of virth between 
stations gives the weight load per station. 

Height of water pressure head. In this problem 
the design water head was the height of the deck 
edge. 

Longitudinal distance between web frames . The 
stresses computed were for the loads on one foot 
of logitudinal distance. Each web frame is assumed 
to support the loads imposed in the lonpitiidinal 
distance between frames. 

Present parameter which represents the number of 
stations from the origin deck centerline. 

Present parameter which represents the number of 
stations from the origin to the deck edge. 

Number of stations from the origin to the bottom 
of the starboard stanchion. 



paJ-i = U 
pa5 = 18 

pa6 = 18 



Number of stations from the origin to the bottom 
of pori stanchion when measuring counter clockwise. 

Number of stations between bottom and top of the 
starboard stanchion. 

Number of stations between bottom and top of the 
port stanchion. 



The following instructions are the sequence control instructions 
and serve to utilize various sub-routines to perform the calculations of 
stresses for a web frame of one particular set of scantlings and load 
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conditions. To alter the web frame and to produce a rational design, 
another set of instructions designated tape number 359-217-222 will be 
required. 



fc 359-217-221 (Name of tape) 

(2k, 6) 

(Selects number arrangement of double length registers. 
This arrangement utilizes 2h binary numbers for the 
mantissa and 6 binary digits for the exponent of the 
number) . 



pi, 


isp gl 


gl forms sin 9, cos 9 and As 


isp 


g2-2 


g2-2 forms horizontal and vertical water loads 

H , V . 
w* w 


isp 


glO-2 


glO-2 forms the total vertical loads by summing 
the values of the weights and vertical water 
loads. 


isp 


g3 


g3 forms sin and cos <f>. 


isp 


gk 


Forms the reciprocal of the moment of inertia 
(l/Ijy jy ), of the section (A^_) and distance 
from the shell to the neutral axis (d^ ^ ). 


OUT 




The following orders are block out instructions 
for transferring information from the case 
memory to the auxiliary drum. 


ca 


c8 


20li.8 (Address of first register of the drum 
storage block). 


si 


9 67 


Selects case output routine to the auxiliary 
drum. 


ca 


cl2 


800 (Number of registers necessary for storage). 


bO 


32 


(Case memory address of the first register being 
blocked, out). 


IN 


isp 


g6-2 


Offsets the x and y coordinates of the neutral 



axis from the shell coordinates. 
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isp 


gl 


Forms 

axis. 


sin 9 and cos 9, and As of neutral 


isp 


g3 


Forms 


sin /, and cos jz i of neutral axis. 


ica 


c6 


+0.0 


+0.0 


its 


tU 




+0.0 


its 


b3 




+0.0 


isp 


g7 


Forms 


m/I at each station due to applied loads 


isp 


g9 


Forms integrals of m^, p^, R and S, into 

a matrix (See Figure VI). 


isp 


g23-2 


Uses C 


1 rout's method to solve the matrix. 


isp 


g2lrl 


Alters 

H *P 0 > 


! V and H to include redundants V-q^, 
R+V and V+S. 


ica 


dl+50 


™0 




its 


tk 




+m 0 


imr 


d3 


1/I 0 


V T o 


its 


b3 




Vo 


isp 


g7 


Forms moments at each station. This moment 
includes the effect of the redundants 



f 



OUT 

ca clO 
si 963 



ca c!7 



bi 332 



ca ell 



2 3k 8 - Drum address of the first register of the 
drum storage. 

Selects "block-in" procedure for transferring data 
from the drum to the case memory. 

+150 - Number of registers being "blocked in." 

Blocks in d , t and t„. 

w s f 

Address of the first register in the core memory 
that will receive the first number moved in from 
the drum. 

+27k8 - Address of the first register of the drum 
storage. 
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si 963 



ca cl8 



bi 732 



IN 

p2, isp g25-l 



Selects "block in" routine from the drum 
storage. 

♦50 - Blocks in 50 registers of information 
of d N,A. 

Address of the first register in the core memoiy 
that will contain d, T . 



Forms C T, <T, , CT and CT for all stations around 
a* b* f s 

the girth., 



The following instructions, called tape 359-217-222, perform the 
stress calculations for a web frame for a given set of scantlings and 
compare the computed stresses with allowable stresses. If the computed 
stresses do not fall within the range of allowable stresses, then the instruc 
tions below will alter the scantlings of the web frame to produce the 
desired stresses. For this program the scantlings are altered by increasing 
or decreasing the width of the flange. Notice that tape 359-217-222 repeats 
all the instructions contained in 359-217-221 and in addition contains a 
method of altering the flange. While this may seem unnecessarily repetitious 
two complete sequence programs were deemed more convenient to handle and less 
susceptible to error than one partial tape with a connective addition. 

fc 359-217-222 

(2U,6) 



Pi, i s P gl 
isp g2-2 
isp glO-2 
isp g3 
isp gli 



gl forms sin 9, cos 9 and As. 
g2-2 forms V and H . 
glO-2 forms 

g3 forms sin / and cos 
gk fonns 1/I N-A- , and 



OUT 



ca 


c8 


20li8 - Drum Address. 




si 


9 67 


Selects drum output routine. 




ca 


cl2 


Ooo - Number of registers being transferred. 


bO 


32 


Core memory, address of first register being 
blocked out. 


IN 






isp 


g6-2 


Forms J„ _ A? ). 




isp 


gl 


sin Q, cos 0, and As of neutral axis. 




isp 


g3 


sin $ and cos <f> or neutral axis, 




ica 


c6 


+0.0 +0.0 




its 


tli 




+0.0 


its 


b3 




+0.0 


isp 


g7 


Forms m /I at each station 
n' 




isp 


g9 


Forms integrals of nr, p . a , R and S. 
Figure VI). *0 O' ^0 


(See 


isp 


g23-2 


Matrix solution. 




isp 


g2l*-l 


Altbrs loads to include redundants. 


p o' 1 


ica 


dl +50 


+n b +m Q 




its 


tti 




+m Q 


imr 


d3 


o 

o 

H 

N 

H 




its 


b3 




V 1 . 


isp 


g7 


Forms moments at each station. 




OUT 






ca 


clO 


23^8 - Drum address. 




si 


963 


Selects "block in" from drum internal routine 


ca c!7 


+15>0 - Number of registers being blocked 


in # 



and 
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bi 332 


Address of first register in core memory to 
receive block in from drum. 


ca ell 


27 l 48 -Drum address of first register in 
,r block in" group. 


si 963 


Selects "block in" internal routine. 


ca cl8 


+50-Number of registers being blocked in. 


bi 732 


Address of first register in core memory to 
receive blocked in data. 


IN 





isp g25 _ l Forms o' , <?, , ^ and <r for all station, 

^ a rr f s 



OUT 




ca cl9 


+2li98-Drum address of first register to be 
blocked in. 


si 963 


Selects "block in" internal routine. 


ca cl8 


+50-Number of registers being blocked in. 


bi I 482 


Core address of first register receiving blocked 
in data. 


aO 976 
cp p2 — 14 


Forms counter '00 limit number of alteration 
cycles to 3(l + number in 976). 


spO 


Stops computer after completion of 3 cycles 



IN 



ica c3 


-1.0 -1.0 


its t3 


- 1.0 


program. 


Note that the above counter was used for the test runs of this 
In any case the alteration procedure should be limited to pre- 



vent the possibility of the computer obtaining an unreasonable answer. As 
an example, should the web frame be to large the web flange can be reduced 
by predetermined amounts. However, large reductions could give a negative 
length of flange. A future programmer may want to change the value used 
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in the counter and therefore the number of cycles e However a counter must 
be used in all cases. 

Now it becomes necessary to compare the values of the computed 
stresses with allowable values. To do this the computed values must be 
plus and a small program must be xjritten to insure that all values are 
plus. Also, the maximum stress value at each point of equal scantlings 
must be determined. Remember that the web frame is symmetrical port and 
starboard, but the values of stress are not symmetrical. After determining 
the maximum value of stress at each station, the values are compared with 
allowable values and the scantlings altered. 



ict pal 
p2, ica e5+c 







icp ml 



if <r is minus perform sub-routine ml. If 
plus, do next instruction. 



its tl 




ica e5>+pb8+pb8 




+ ^(port side). 



icp ml 



Checks sign of CT, a 



its t2 



icp m2 



isu tl 



-^(starboard) C^(port) - CT"( starboard) 

Checks sign of (7^. (port) - 0^ { starboard) and if 
minus, goes to m2. 



ica t2 





p3, isu c20 



-(+20,000) CT -(+20,000) (Lower value of stress 



limit) . 



If ((T-20,000) is minus, go to m3 to reduce 
scantling values. 



icp m3 
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isu c21 

icp pit 

ica b3+c 

iad cl 

its b3+c 
ics c3 
its t3 
pU, OUT 
ca p2+3 
su c7 
ts p2+3 
IN 

ict p2 
ica t3 
icp p6 
OUT 
ca cl9 
si 967 
ca cl8 
bo I482 

ca c8 

si 963 
ca cl2 



-8000 (CT-20 o 000) -8000 (8000 /in 2 

is the range of stresses). 

If value is minus, perform sub-routine pl| which 
moves forward to the next station. 



+w f 


+W f 




w +.5 (Increases flange width). 




w +,5 




f 


-(- 1 . 0 ) 


+1.0 




+1.0 



ica e5>+pb8+pb8 ica eE>+pb8+pb8 

-2 ica e5>+pb8+pb8-2 

ica e£+pb8+pb8-2 



- 1,0 



- 1,0 



+2lt98-Drum address register. 

Selects drum block out procedure. 

+50-Number of registers being blocked out. 

Core memory address of first register of group 
being blocked out, 

+20li8-Drum address register of first number in 
a group being blocked in. 

Selects drum block in procedure, 

800-Number of registers being blocked in. 
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ca cl2 


+800 


+800 


ad cl7 


+150 


+930 


ad c2 


+30 


+980 


ts t8 






ca t8 


+9 80 




si 963 


Selects 


drum block in procedure 


ca cl2 


+800 




bi 980 







The two block in procedures put numbers and instructions 
back into the core memory in preparation for another cycle of stress 
calculations of the altered scantlings® The block out instruction 
prior to the block in instructions puts new values of web flange 
width into the columns of numbers on the drum. The first block in 
instruction moves numbers into the core memory® The second block 
in instruction came as a result of "debugging" operations when it 
was determined necessary to replace all the addresses in the in- 
structions back to their starting value involving the port side it 
was necessary for the programmer to write instructions that would 
decrease the addresses of the registers containing the starboard side 
values. In other words it was necessary to move up a column of 
numbers in the core memory. At the beginning of the second cycle 
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the instruction addresses were in error by a number equal to the number 
of cycles in that particular sub-routine. By blocking in the original 
instructions which remain on the drum into their assigned positions 
all the addresses were reset to their original value. 

For cycles within a basic cycle of operations it is necessary 
to write a reset routine that returns the addresses back to their 
original values at the end of each inner cycle. 



IN 

isp gif. 
isp pl+11 



p6, i STOP 



Starts cycle at gU. 

These two instructions skip the block out 
procedure in the main cycle. Since the block 
out has occured during the first cycle, the 
values remain in the drum and need not be 
blocked out again. 



The following sub-routines are used in conjunction with 
the comparison and alteration sequence. 
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ml, ita ml+2 



imr c3 
isp 


x-loO 


(-l e 0x-(^) = + (7^ 
(Makes CT^ plus). 




m2, ica tl 


♦<£ 


+0^( starboard) 




isp p3 


m3, ica b3+c 


+w f 


+w 




isu cl 


~=5 


W f “o5 




its b3+c 




w -.3 


isp p3+7 


d9. 


+2lr98-Drura address. 




c20. 


+20,000 . 0-Lower 
value. 


limit of allowable 


stress 


c21. 


+8, 000 .0-Range 


of allowable stress 


values. 



Note that for this problem a very wide range of values was 
used to prevent the possibility of the computer seeking 
an Inpossible solution* If the stress change caused by 
the flange alteration were to be larger than the stress 
range, then a solution would be impossible with this program. 



976 -2 

i START AT pi 

In programming it is necessary to use various constants to change 
units, to form equations and to change generated values by arithmetic 
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manipulations. These constants must be assigned registers or flads to 
be of use in the program. The constants listed below were assigned flad 
cl, c2, etc., the c being useful in remembering the nature of the input. 

Tape fc 359-217-101 
(25,6) 



cl. 


+0.0 






c2. 


+ 30.0 


A constant used to form the width 
plating which is considered to be 
forming the web beam. 


of shell 
effective in 


c3. 


-1.0 


A constant useful in changing the 
value in the MIA. 


sign of a 


c5. 


+1.5 


Used in calculation of the moment 


of inertia. 


c5. 


+65.0 


2 

Weight of a cubic foot of sea water in lbs/ft . 


c6. 


+0.0 


Used to clear registers where the 
register must be zero. 


value in the 


c7. 


+2 


Whirlwind number used in changing 
interpreted mode. 


addresses in the 


c8, 

c9, 

clO, 

ell, 


+2058 

2850 

+2358 

+2758 


A Whirlwind number used to assign 
the auxiliary drum. 


addresses on 


cl2, 

cl3. 


+800 

+216 


Number of registers in a blockout 
procedure. 




oH>» 


+1 


Whirlwind number used in changing 


addresses. 


clS, 


+12.0 


Converts feet or inches to inches 


or feet. 


cl6. 


+0.707 


1/ 2 = sin 55° or cos 55°o 




cl7, 
cl 8, 


+i5o 

+5o 


Number of registers in block out 
procedures- 





tl, through t8, +0.0 - Temporary registers for calculations. 



92 



The following preset parameters are part of the inpu-.s t' :.t 
do not vary with vessel type . Note that each preset parameter lists' below 
is generated from the present parameters in the variable input tape. 



pbl = pal-1 


Number of segments in the half girth of a 
web frame. 


pb2 = pal-2 


Tiro less than the number of stations on a half girth 
Used in moment formation. 


pb3 = pa2-l 


Number of segments between the origin and deck edge. 


pbli = pa3 +pa3 -2 


Number of registers from the origin register to 
the bottom of the starboard stanchion. 



Also contained in Tape 359-217-101 are the storage assign- 
ments for generated numbers. 



132 1 al, +0.0 
DITTO 

U8 +o.o 


Sine of the angle storage c Note that this 
instruction puts zeros in 25 double length 
registers* 


182 | a2, +0.0 
DITTO 
58 +0.0 


Cosine storage e 


232 | a5. +0.0 

DITTO 
58 +0.0 


Girth segment lengths storage 


532 | eU, + 0,0 
DITTO 
58 +0.0 


Storage for vertical water loads. 


582 | e3, +0.0 
DITTO 
58 +0.0 


Storage for horizontal water loads < 


832 | d3, +0.0 
DITTO 
58 +0,0 


Storage for the inverse of the moment of inertia 
1 f l 

1 X N.A, 























! I 
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882 | d2, +0.0 
DITTO 
1*8 + 0.0 

7321 dl, +0.0 
DITTO 
1*8 + 0.0 

?82| si, Library'’ sub-routine for finding square roots. 

The remainder of tape 339-217-101 consists of the various 
sub-routines that perform the functions listed in the flow diagram. 

Figure V„ A discussion of each sub- routine is placed at the beginning of 
the list of instructions that comprise the sub-routine® 

The form of presentation of the follotiring sub-routines was 
found to be an aid in the organization and understanding of the problem. 
Since continuity of thought is so important to the writing of a work- 
able program, the increased effort in writing the flow of information 
into the MRA, accumulation of data in the KRA and the values trans- 
ferred from the MRA was considered worthwhile. Other methods of 
organizing material could have been used; however, the authors found 
that this method was both effective and presented the ideas in a 
convenient form. 
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The foil wing sub-routine forms the sin 0 and cos 9 of the 
frame segments by the use of the basic equations listed below. 



As 


= Ax^ + 


Ay 2 


n 


n 




Ax 


= x -x 




n 


n n- 


-1 


% 


= YTT 
n n- 


-1 


sin 0 


= Ay/As 




cos 0 


= Ax/As 






Instruction 


Input to MRA 


In MRA Transfer 








from MRA 


gl. 


ita gl+21 


This instruction transfers the address in the 






control element to the 


address of the instruc- 






tion in gl+21 and adds 
instructions) , 


one. (See glossary of 


icr 


pbl 


Sets the cycle count to the number of segments 
on the side of the frame. This routine will be 






performed pbl times 




ica xl+2+c 


♦x 


+x 






n 


n 


isu 


xl+c 




x -x n 






n-1 


n n-1 


its 


tl 




Ax 

n 


imr 


tl 


X Ax 


Ax 2 






n 


n 


its 


t2 




Ax 2 



n 



9S 



Instruction Input to MRA In HR A Transfer 

from MRA 









(continued) 


ica 


yl+2+c 


+y n 


+y n 


isu 


yl+c 


" y n-l 


VVi 


its 


t3 






imr 


t3 


<] 


Ay 2 


iad 


t2 


2 

+Ax 

n 


Ay^ + Ax' 
•'n i 



SP si si is a sub-routine that extracts the square root 
of the number in the MRA and returns the square 
root of the number back into MRA. 



its 


a5+c 




As = v/Ax 


+ A / 


ica 


tl 


+Ax 

n 


+Ax 

n 




idv 


a5+c 


-Ax n 


Ax /As 
n n 




its 


a2+2+c 




(cos 0) 


- Ax /As 
n n 


ica 


t3 


+Ay n 


O 




idv 


a5+c 


vAs 

n 


Ay /As 




its 


al+2-*-c 




(sin 9) 

n 


= Ay /As 
J n' n 


ict 


gl+2 


Puts gl+2 in the 


control element and starts the 





cycle over again. 



isp 



Address put into this instruction by the ita 
instruction in flad gl. 
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The sub-routine that follows forms the water loads from the 
origin to the main deck a In this sub-routine the water loads are assumed 
to be uniformly distributed along the shell segment and the center of 
pressure is assumed to act at the mid-point of the segment. After formation 
of the total load on the segment, the load is divided by 2 and half of the 
load is placed at each end of the segment. The approximation of a con- 
tinuous beam under a trapezoidal shaped pressure load by a uniform load 
is only valid when the shell segments are small in size. As the loads 
are formed, the program converts the loads into horizontal and vertical 
components. 

z 




Instruction Input to MRA In HRA Transfer 

from KRA 

ita g2+l8 (See glossary of instructions). 

icr pb3 Sets cycle counter for number of segments in the 
starboard side from origin to deck edge. 

g2, ics yl+c -y n -y n 

n+1 -W y „ " 



isu yl+2+c -y 
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Instruction 



Input to MRA 



In HRA 



Transfer 
fron HRA 



(continued) 



imr 


cl 


o.5 


- (y n.lV /2 * - y av 


iad 


e2 


+H(water head) (ft) 


H - y av 


imr 


c5 


a 6k ./ft 3 




imr 


a5+c 


*As (ft) 


6hAs(H-jJ . F w 


imr 


cl 


*0.5 


F ^/2 (assumes simple support) 


its 


tl 




V 2 


imr 


a2+2+c 


* (cos 9) -a 
n+1 


(cos eJ^F/2 . (V w ) n+1 


its 


eU+2+c 




( Vn*l 


iad 


eii+c 


V 

w n-l 


V +V « V 

w . n w 

n-1 


its 


elj+c 




<Vn 


ica 


tl 




* ( V 2 >n 

F 


imr 


al+2+c 


+(sin 9) 

n 


( w sin 0) = (H ) 

c. n n 


its 


e3+2+c 




(H ) 
w n 


iad 


e3+c 


(H ) 
w n-1 


(H ) _ + (H ) . H 

w n-1 w n w 

n 


its 


e3+c 




(H ) 
w n 


ict 


g2 


cycle back to g2. 




isp 




Address provided by ita instruction. 



98 , 

Formation of vertical forces is accomplished by the programs 
listed below. Vertical forces are the algebraic addition of buoyant forces 
which act in the upward direction and the weight forces which act down- 
ward. Upward forces are taken as positive. Note that a separate program 

is required for each side. The V forces are symmetrical and can be 

w 

stored in !?0 registers. The weight forces are non- symme trical and the 

weight for each of the stations must be stored. In the program for 

forming the starboard vertical forces, the addresses of the V moves down 

w 

a column of registers of increasing addresses which is the manner of 
operation of the eli+c address. To move the other direction up a column 
of numbers a small Whirlwind program is required to reduce the addresses 
by 2 for each cycle. 

Instruction Input to HRA In HRA Transfer 

from MRA 



ita 


(See glossary of instructions) 


• 


icr 


pa2 






gio. 


ica eli+c 


+V 


+V 




w 


w 


iad 


e5+c 




+v -vr = v 

w t t 


its 


e£+c 


End of starboard 


side program. 


ica 


eU+c 


+V 


+V 




w 


W 


iad 


e5+pb8+pb8 


- w t 


+V -W x = V 
w t t 


its 


e5+pb8+pb8 







Goes out into Whirlwind to change addresses 
of weight loads. 



OUT 
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Instruction 



ca glO+ii 
su c? 
ts glO+li 
td glO+5 
IN 



ict glO+2 
isp 



Input to KRA In KRA Transfer 

from MIA. 



(continued) 

iad e3>+pb8+pb8 iad e5>+pb8+pb8 
-2 iad e5+pb8+pb8-2 

iad e5>+pb8+pb8-2 

e5+pb8+pb8-2 



After forming sin 0 and cos 0 of each of the segments, it was 
necessary to find the sin and cos of the angle / formed by the tangent to 
the curve at each station and the horizontal axis*, cos j£ and sin $ were 
determined by averaging the cosines and sines of the angles of the chords 
on either side of the station. While it was considered to be more accurate 
to find the tan or the slope of the curve at each station, the possibility 
of tan $ = Ay/Ax becoming infinity rendered this impossible. Further, only 
sin / and cos / were necessary for the purposes of the program and tan ji 
would have required additional storage and programming. This implies two 
conditions that do not exist. One p the ’"ate of change of angles should 
be equal and secondly, the 



sin 



(e 1 + © 2 ) 



does not equal 



sin 0^ + sin © 2 



2 



does not equal 



cos 9^ + cos 0£ 

2 * 

However, the resulting error will he small as long as the station points 
are judiciously selected to maintain small nearly equal angle changes. 
Should more accurate angles be necessary the programmer could utilize 
trigonometric identities and equal angle changes. 

An additional problem in this sub-routine was to insure that 
the centerline values of angles on deck and bottom represented 0° and 180° 
respectively. As was mentioned earlier the sin and cos of 135>° were placed 
at the corner station. 



Instruction 


Input to MRA 


In MRA Transfer 

f^om MRA 


S3, 


ita g3+20 


(See glossary 


of terms ) . 


icr 


pb2 






ica 


al+2+c 


♦sin 9^ 


+sin 0^ 


iad 


al+lj.+c 


♦sin 0£ 


+sin 9^ + sin 9^ 


imr 


cl 


xo.5 


sin 9 + sin 9 0 

sin i 


its 


al+2+c 




sin <f> 


ica 


a2+2+c 


+cos 9^ 


+cos 9^ 


iad 


a2+li+c 


+cos 0 2 


cos 9^ + cos 9^ 
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Instruction Input to MRA. In MRA Transfer 

from MRA 

(continued) 

cos 0^ + cos ©2 

imr cl. 0.5 5 = cos <f> 



its 

ict 



ica 

its 

its 

ics 

its 

ica 

its 

ica 

its 

ics 

its 

isp 



a2+2+c cos <f> 

g>2 

The following orders sets 0° and l80° in the 
center line deck and keel and, also sets 135 
in the corner angle. 

c6 +0.0 +0.0 

al +0.0 



al+pbl+pbl +0.0 

c3 +1.0 +1.0 

a2 +1.0 



c3 -1.0 

a2+pbl+pbl 

cl6 +0.707 

al+pb3+pb3 

cl6 -0.707 

a2+pb3+pb3 



- 1.0 



+0.707 



-0.707 



- 1.0 



+0.707 



-0.707 



Since web frames are built up sections it is necessary to 
compute the moment of inertia at each station around the girth. The width 
of shell plating assumed to be effective in forming the web beam is 
thirty times the plate thickness. 
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The following sub-routine uses the relationships below in 
forming the area, moment of inertia and distance to the neutral axis 
of the web from the shell* 



A r ■ Vf 



m f - (d » * 2“ )A r 



- (d „ * “ )m f 



A = d t 
w w w 



m = ■ A 

w 2 w 



L t d 

r w w w 

w = 1^5 m w = ~ 



A = 30 t 
s s 



m = —A 
s 2 s 



t 

T S 

1 = -T—p m 

s 1.5 s 



A, = A +A +A 
t f w s 



m, = m„+m +m 
t f w s 



VI +1 

f w s 



m, 



-A t 



NA 



I t" A t d WA 
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Instruction 


Input to MRA 


In MRA 


Transfer 
from MRA 


gU, 


ita g$+2h 








icr 


pal 








ica 


bl+c 




+d 

w 




irar 


b2+c 




d x t - A 

W V V 




its 


tl 






A 

V 


imr 


bl+c 




A v*4* 




irar 


cl 


*+o.$ 


0*5 * \ * d w 


» m 
w 


its 


t2 






ra 

w 


ica 


bl|+c 


+t f 


+t f 




imr 


b3+c 


X 


tf * W f " A f 




its 


t3 






A r 


ica 


bli+c 




+t f 




imr 


cl 


xo.5 


0.5 A t^ 




iad 


bl+c 


♦d 

W 


V 2 




imr 


t3 


XA f 


A f (d „ * V 2) 


- j* f 


its 








ra f 


ica 


b$+c 


+t 

s 


+t a 




imr 


b5+c 


*t 

s 


<v 2 




imr 


c2 


X30 


30(t ) 2 . A 

S 8 




its 


t5 






A 

Q 



loU 



Instruction 


Input to MRA 


In MRA Transfer 

from MRA 






(continued) 


imr 


b5+c 




A A t 
s s 


imr 


cl 




A t 
s s 

2 = m s 


its 


t6 




m 

s 


ica 


t2 


+m 

>7 


m 

w 


iad 


t h 




m + m„ 
w f 


isu 


t6 


-in 

s 


m +m -m * M 
w f s t 


its 


t7 




"t 


ica 


tl 


+A 

w 


+A 

w 


iad 


t3 


+A f 


A +A 
w f 


iad 




+A 

s 


A +A +A jn xr A. 
s w f t 


its 


d2+c 




A t 


ica 


t7 


♦ M t 


* M t 


idv 


d2+c 


t* 


M t 

a; - Va. 


its 


dl+c 




d M.A. 


g5, ; 


ica t2 


+m 

w 


w 


imr 


bl+c 


Xd 

w 


Ad 2 

m d = 
w w 2 


idv 


cl; 


1.5 


A d 2 

W W -j- 

0 **• JL 

3 w 



its tl 
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Instruction 

ica bli+c 
imr cl 
iad bl+c 
imr tb 
iad tl 
its tl 
ica t6 

idv ci-j. 

imr b5+c 

iad tl 

its tl 

ics d2+c 

imr dl+c 

imr dl+c 

iad tl 

its tl 

ics c3 
idv tl 



Input to MRA 


In MRA 


Transfer 
from MRA 


(Continued) 










*o«5 


t f /2 




d 

w 


W 2 




X rn f 


m f ( d w +.5t f ) 


- VW 2) 


+1 

w 


I +I ^ 
w f 


( W 


+m 

s 


+m 

s 




-ri.5 


,2 

m w t 

s s s 

1.5 “ 3 




x t 

s 


3 

w t 

S S t 

** J- 

3 s 






W J s - 


^about shell 






I shell 


~ A t 






Kd N c A. 


" AI d N.A. 




* d N.A. 


- AI 4.A. 




^about shell 


2 

^.A. = I shell“ AI ^I.A. 








-(-1.0) 


+lo0 




fI N 0 A. 


1//:C N.A t 








Its 



d3+c 
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When working with beams, all forces and moments are considered 
to act along the neutral axis. To make the program more general, the 
original coordinate inputs were the x and y coordinates of the inner 
surface of hull shell. The coordinates of the neutral axis may be formed 
by offsets from the shell coordinates. The reader will notice that such 
a program allows alteration of the web section as desired,, to produce the 
practical results. The offset distance is ^ 




^.A. * x s ' Va. Sln 
y H.A. ' y s * d N.A, 005 

The following program forms the coordinates of the neutral axis. 
Since this program is a continuation of the above sub-routine, the counter 
instructions are still in effect. 



Instruction 



Input to MRA 



g6, ics dl+c +d^ ^ 

idv cl5 4-12 



In MRA. 



Transfer 
from MRA 



- d W.A. (il,ch) 
-d N _ A A2 (ft) 



"%.A.^ 



its tl 
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Instruction 


Input to MRA 


In MIA. 


Transfer 
from MRA 






(continued) 




imr 


al+c 


* siJrl ''shell 


d H.A. sin * 
12 




iad 


xl+c 


+ x shell 


\ 

5 

CO 

0 

<3 

0 

"d 


s 12 


" N.A. 


its 


xl+c 






X N.A. 


ics 


tl 


CM 

< 

o 

< 

0 

1 

i 


+d N e A/ 12 




imr 


a2+c 


XC0S '‘shell 


cos i 




iad 


yl+c 


7 shell 


d„ . 

N 0 A„ , 

cos * - J'n.a. 


its 


yl+c 






7 n.a. 


ict 


gk 








isp 











The moments at each station are formed by multiplying the 
summation of the horizontal and vertical components of the loads by the 
vertical and horizontal projections of the segment* 
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(x -x , ) 
n n-1 




( Wi 



) + m 

n 



0 



The following sub-routine forms the moments at each station due 
to the weight and water pressure loads on the starboard side. Note that 
the first few instructions put +0.0 into tl, t2, t3. This was necessary 
due to their use in forming the summations, In addition when using this 
sub-routine to form m^, tlj. and b3 must be set to zero. When using this to 
compute the bending moment after solving for the redundant, the value of 
the redundant must be put into tii and m^/l^ in b3. 



Instruction 



Input to MRA 



In MRA 



Transfer 
from MRA 



g7, ita g8+35 

icr pbl 

ica c6 

its tl 

its t2 

its t3 

ica e£+c 

iad tl 



+ 0.0 



hV T 

n-1 
n~2 

+ > V 



+ 0.0 



Vi 

n-1 



V 



+ 0.0 

+ 0.0 

+ 0.0 



its tl 



> 



n-1 

—7 . 



ics xl+2+c 



-(x ) 

na n 



(-x ) 

na n 



imr tl 



n-1 

“7 



V 



n-1 
Ax > V 



Z’ 



(m ) 
n v 
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Ins traction Input to MRA In MRA Transfer 

from MRA 

(continued) 



its 


t3 






(m ) 
n 


ica 


e3+c 


♦H 

VI 


+H 

w 




iad 


t2 


n-2 

y« 

Z_ ^ w 


n-l 

Z ^ w 


n-l 


its 


t2 






\ ' H 

/ s w 


ics 


y+2+c 


^” y N.A. 


^“ y N.A. 




iad 


yl+c 


^N.A^n-l 


~ f y N.A,^N “ 


(y N.A.Vl] 


imr 


t2 


n-1 

v \ H 
/ \ w 


n-l 

A 'V 7 H 
J w 




iad 


t3 


(m ) 1 

n n-1 


(ra ) 
n n 




iad 


t b 


n-l 

X<"n> 


(m ) 
n 




its 


tii 






( ”n ) 


imr 


d3+2+c 


X 1/1 


m /I 
n 




its 


b3+2+c 






V 1 


ict 


g7+8 









no. 



The following section of the sub-routine forms moments port 
side. Remember that the coordinates of neutral axis and the horizontal 
loads are stored for one side. The sub-routine must move up the 
stored column of numbers to form the moments. 



Instruction 


Input to MRA 


In MRA 


Transfer 
from MBA 


icr pal 








g8, ica e5+pal+pal+c ♦V 


♦V 

n 






n-2 


n-1 


n-1 


iad t3 


X* 


Z V ‘ V n 


-X' 








n-1 


its tl 






x* 


ica xl+pb2+pb2 


+X n+1 


♦X 

n+1 




isu xl+pbl+pbl 


~ x n 


n+1 n 


* A 

X 




n-1 


n-1 




imr tl 


X v 


4 * 2 t 


» TTl 

n 


its t3 






(m ) 








n 


ica e3+pbl+pbl 


♦H _ 
n-1 


♦H _ 
n-1 






n-2 


n-2 


n-1 


iad t2 


2-- 


H +H 'u 

/ x n-2 n-1 / v n 



HI. 



Instruction 


Input to MRA 


In MRA 


Transfer 
Troiri KRA 






(continued) 




its 


t2 






n-1 


ica 


yl+pb2+pb2 


^ +y N.A.^n+l 


^ +y N,A. ^n+1 




isu 


yl+pbl+pbl 




^W.A^n+l 


“ ^ 7 K,Ajn] 






n-1 


n-1 

-s / 





imr t2 



"l 



H = (m n ) H 



iad t3 



(+m ) 

n v 



(m ) + (m) = m 

n v n n n 



iad ml+pbl+pbl+c 



n-1 
7 



m 



n 



n-1 



n 



m + m = 
n 



m 



imr d3+pbl+pbl l/I 



n 



m/I 



OUT The program now goes out into Whirlwind to 

change the addresses for each cycle. 



ca g8+3 
td g8+U 



ica xl+pb2+pb2 



ica xl+pb2+pb2 



Puts the address of gS+3 into the address of 
g8+lj now reads isu xl+pb2+pb2 , 



su c7 
ts g8+3 
ca g8+7 



-(+2W.WJ 



ica e3+pbl+pbl 



ica xl '-pb2+pb2-2 

ica xl+pb2+pb2-2 



ica e3+pbl+pbl 
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Instruction Input to MM In MRA Transfer 

from MRA 

(continued) 



su 


c7 


-(+2 ¥.¥.) 


ica e3+pbl+pbl-2 


ts 


g8+7 




ica e3+pbl+pbl-2 


ca 


g8+10 


ica yl+pb2+pb2 


ica yl+pb2+pb2 


td 


g8+ll 


(Instruction in g8+ll 


is now isu yl+pb2+pb2). 


su 


c7 


-(+2 ¥.¥.) 


ica yl+pb2+pb2-2 


ts 


g8+10 




ica yl+pb2+pb2-2 


ca 


g8+l5 


+imr d3+pbl+pbl 


+imr d3+pbl+pbl 


su 


c7 


-(+2 ¥.¥.) 


+imr d3+pbl+pbl-2 


ts 


g8+l5 




imr d3+pbl+pbl-2 


IN 






ict 


g8 






isp 









Since all moments due to applied loads and forces had been 
formed, it was possible to perform the integrations of moments to obtain 
the slopes and deflections previously discussed in the procedure. At 
this point in the programming the redundants m^p^, and must be con- 
sidered. m^ will appear as a constant value of moment at each station 
throughout the girth. The moments due to p^ and depend upon the 
value of y and x at each station. Since the web frame was assumed to 

m , , , = p v- + q x + m 

redundant 0 n 0 n 0 
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be symmetrical* the slope changes and deflections around the entire girth 
caused by the redundants will be double the slopes and deflections on one 
side, Fo: convenience in integration, the redundants will be integrated 
around the starboard side. 

The method of integration used was the trapezoidal rule which 
implies linear slope of the moment curve between stations on the moment 
curve. The area under the curve 




n=l 



(H) ] 

l iViJ 




A 




(H) 

'I'2 







As can be seen from the figure, the trapezoidal rule will cause small 
error providing the curve is smooth and stations are closely spaced. 
Despite the small error, the above method of integration was chosen 
because it is readily applicable to computer programming and because 
random station spacing was necessary to approximate a curved figure 
by a series of straight lines. 



In the process of "debugging" the program several improvements 
in the sequence of operations for the formation of the integrals were 
made. The results of these improvements and errors in addresses of the 
initial program for the integration procedures alters the procedure as 
written. The techniques used and explained in the following pages are 
correct. 

The major improvement was the inclusion of the integration 
routines under one sub-routine. All the integrals for iHq, p^, and 
m /I were performed under one cycle counter. Each of the stanchion 
redundant forces R and S were under a separate cycle. The routine 




was not performed. To obtain the coefficient of p^, the coefficient of 
m Q which is 

was used. The group of instructions included under flad g 9 were added 
to insure that zeros were initially placed in the registers that were 
to hold the matrix. Also, the addresses for forming the port side 
integrals were found to be in error by two registers. 

The instructions listed below are correct and replace the 
integration routines listed under flads gll through g22 of the detailed 
description. Unf ortunately, these alterations were not made in sufficient 
time to allow the correction of the detailed explanation. These instruc- 
tions should be read from left to right. 



n5 



g9, ita g22 


icr 30 


ica c6 


its dl+c 


ict g2+2 


icr pbl 


gll, ica d3+c 


iad d3+2+c 


±rar a5+c 


iad dl 


its dl 


gl2, icaxl+c 


imr xl+c 


imr d3+c 


its tl 


ica rl+2+c 


imr xl+2+c 


imr d3+2+c 


iad tl 


imr ap+c 


iad dl+2l; 


its dl+2li 


gl3, ica yl+c 


imr d3+c 


its tl 


ica yl+2+c 


imr d3+2+c 


its t2 


iad tl 


imr a5>+c 


iad dl+2 


its dl+2 


glli, ica tl 


imr yl+c 


its tl 


ica t2 


imr yl+2+c 


iad tl 


imr a3+c 


iad dl+12 


its dl+12 


gl5, ica b3+c 


iad b3+2+c 


imr a5+c 


imr cl 


iad dl+5>0 


its dl+^O 


gl6, ica yl+c 


imr b3+c 


its t3 


ica yl+2+c 


imr b3+2+c 


iad t3 


imr a5>+c 


imr cl 


iad dl+52 


its dl+^2 


gl7* ica xl+c 


imr b3+c 


its t3 


ica xl+2+c 


imr b3+2+c 


iad t3 


imr a5+c 


imr cl 


iad dl+51i 


its dl+51* 


gl8, icab3+pbl+pbl+c 


iad b3+pal+pal+c 


imr a5+pb2+pb2 


imr cl 


iad dl+5>0 


its dl+50 


gl9, icayl+pbl+pblimrb3+pbl+pbl+c 


its tl 


icayl+pb2+pb2 


imrb 3+pal+pal+ c 


iad tl 


imra3+pb2+pb2 


imr cl 


iad dl+5>2 


its dl+52 


ic sxl+pbl+pbl 


imrb 3+pbl+pbl+ c 


its tl 


icsxi+pb2+pb2 


imrb 3+pal-f pal *• c 


iad tl 


imra!?+pb2+pb2 


imr cl 


iad dl+5l 


its <n.+5h 


OUT 


ca gl8+2 


su c7 


ts gl8+2 


ts gl9+6 


ts gl9+l6 


ca gl9+3 


ts gl9 


su c7 


ts gl9+3 


ca gl9+13 


ts gl9+10 


su c7 


ts gl9+13 


IN 


ict gll 


ica dl+2 


its dl+10 


g20,icr pa5 


icaxl+pbli+c 


isu xl+pbU 


its t3 


irira3+pbl|+c 


imr cl 


its tl 


imrb3+pbL.+c 


iad dl+56 


its dl+56 


icaxl+pbIi+2+c 


isuxl+pbij 


its tli 


imra5+pbli+c 


imr cl 


its t2 


imrb 3+pblt+ 2+ c 


iad dl+36 


its dl+^6 


ica tl 


imrd3+pbli+c 


its tl 


imiyl+pbli+c 


iad dl +16 


its dl +16 
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ica t2 


imrd3+pbli+2+c 


its t2 


imiyl+pbl;+2+c 


iad dl+16 


its dl+16 


ica tl 


imrjl+pbl;+c 


lad dl+26 


its dl+26 


ica t2 


imrxl+pbl;+2+c 


iad dl+26 


its dl+26 


ica t3 


imr tl 


iad dl+36 


its dl+36 


ica tU 


imr t2 


iad dl+36 


its dl+36 


ica tl 


iad t2 


iad dl+6 


its dl+6 


ict g20+l 


ica dl+6 


its dl+30 


ica dl+l6 


its dl+32 


ica dl+26 


its dl+3li 


g21, icr pa6 


icaxl+pb5+c 


isuxl+pbl? 


its t3 


imra5>+pb!?+c 


imr cl 


its tl 


imrb 3+pb6+pb6 


iad dl+58 


its dl+58 


icaxl+pb3>+2+c 


isuxl+pb5 


its tl; 


imra5+pb5+c 


imr cl 


its t2 


imrb 3+pb6+pb6-2 


iad dl+38 


its dl+58 


OUT 


ca g21+l6 


ts g21+7 


su c7 


ts g21+l6 


IN 


ica tl 


imrd3+pb5+c 


its tl 


imiyl+pb5+c 


iad dl+18 


its dl+18 


ica t2 


imrd3+pb5>+2+c 


its t2 


imryl+pb 5+ 2+ c 


iad dl+18 


its dl+18 


ics tl 


imrxl+pb5+c 


iad dl+28 


its dl+28 


ics t2 


imrxl+pb3+2+c 


iad dl+28 


its dl+28 


ica t3 


imr tl 


iad dl+l;8 


its dl+I;8 


ica tl; 


imr t2 


iad dl+l|8 


its dl+Ii8 


ica tl 


iad t2 


iad dl+8 


its dl+8 


ict g21+l 


ica dl+8 


its dl+liO 


ica dl+18 


itx dl+l;2 


ica dl+28 


its dl+l|li 


g22, ispO 
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The results of each integration are stored to form the matrix 
as shown by Fig. VI, The following sub-routine performs the integration 
of the change of slope around the girth of the web frame due to redundant 
moment m^. 



A0 



■/ -r ds - 2 f t ds - 2 n o/ ( f * rt ) ^ • m o/ ( y- ♦ rt )ds „ 
o o o n n+1 o n n+1 



Instruction 



Input to MRA 



In MRA 



Transfer 
from MRA 



gll, ita gl5+9 (See glossary of instructions). 

icr pbl (Sets cycle counter equal to the number of - 

stations on a side minus l). 



ica 


d3+c 


+1/1 . 
' n-1 


+1/1 
' n- 


•1 


iad 


d3+2+c 


+1/1 

n 


+1/1 
' n 


+ 1/1 . 
' n-1 


imr 


a5+c 


xds n • 


(1/1 
' n 


+ l/I _ 
' n-1 


iad 


dl 


n-1 

*Z (1/I n * 


l/I )ds n 
' n-1' n-1 y 


n 

> ' 



its dl 



n 

Z (l/I + 1/1 , )ds 

' n n-1 n 



The instructions below find the deflection in the y-direction 
due to the moment caused by the shear redundant q^. 



0 A n-1 

Ay = J' x(q () x)ds = 2 J x(q Q x)ds « 2q 0 ^~^ 
0 0 



(x n) 2 (x ) 2 
n-1 n 



n-1 



n 



ds 

2 
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Instruction 


Input to MRA 


In MRA Transfer 

from MRA 


gl2. 


ica xl+c 


+ Vi 


+x n 
n-1 


irar 


xl+c 


* X n-1 


<Vi> 2 


imr 


d>c 


* 1 / I n-1 


^n-l'* / I n-1 


its 


tl 




(x -,) 2 /l n 
n-1 ' n-1 


ica 


xl+2+c 


+x 

n 


X 

n 


imr 


xl+2+c 


11 *n 


2 

X 

n 


imr 


d3+2+c 


xl/^ 


(x) 2 / I 

n n 


iad 


tl 




(x ) 2 /I + (x n ) 2 /I - 
n ' n n-1 ' n-1 


imr 


a5+c 


x ds n 


(x ) 2 /I + (x _ ) 2 / I _ ds 

n ' n n-1 ' n-1 n 






n-l 


n 


iad 


dl+2l| 


+ 'N ds , (x /I ) 2 
/ v n-1 n' n 


+ (x , ) 2 /I _ \ (x ) 2 /I +(x 
n-1 ' n-1 / s n ' n n- 


its 


dl+2li 




n 

(x) 2 /l+(x -)' 

£ x n n n-1 



The deflection in the x-direction due to m^ was found by per- 



forming the integration 



Ax 



0 A n-1 

/ y ds. 2 J |(y n -l/ I n-l K( 5'n /I n ) 

0 0 1 



~i ds 



n 



The instructions starting with gl3, perform this integration, 
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Instruction 


Input to MRA 


In MRA Transfer 

from MRA. 


gl3. 


ica yl+c 


ty n-l 


%.i 


jjnr 


d3+c 


x ^n-l 


VXn-l 


its 


tl 




y n-l^n-l 


ica 


yl+2+c 


+y n 


+y n 


imr 


d3+2+c 


X i/i 
7 n 


hA. 


iad 


tl 


* -^n-l^n-l 


V.'vAi’ 


imr 


a5+c 


X ds n 


(y /I +y /I )ds 

n n ■ 7 n-l / n-1 n 






n-1 


n 



+ / (y /I +y -i/I )ds _ / (y /I 4-y Ar )ds 

X v rv n n-1' n-1 n-1 X . n' n ^ n-1' n-1 n 



1 1 



its dl+2 



n 




The following section of the sub-routine performs the integration 
of the moment due to = p^ around the girth. 



AS 






n+1 

n+1 



ds 

n 

2 
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Instruction Input to MRA In MRA Transfer 

from MRA 



gilt. 


ica yl+c 


^n-1 


+y n-l 


imr 


d3+c 


X 1 / I n-1 


y n-l/ I n-l 


its 


tl 




y’n-l^n-l 


ica 


yl+2+c 


+y n 


+y n 


imr 


d3+2+c 


X!/I n 


y n /J n 


its 


t2 




y /i 

J n / n 


iad 


tl 


+y n-l //l n-l 


y /I +y _/l _ 

J n n J n-1' n-1 


imr 


a3+c 


X ds n 


(y /I +y ../I , )ds 
n n n-1 n-1 n 






n-1 


n 


iad 


dl+10 


5 y n-l / ' I n-l +y n- 


-2 //l n-2^ dS n-l / .^n^V^n-l^n-l 



its dLrt0 zlwv/n-l^n 



The x-deflection due to the axial redundant force was 





ds = 2p 0 




+ 



( y n -i r 




ds 

n 

2 



The instructions given in gl£, performs the p^ moment x-deflection. Note 
that the values of y ^/l n _^ y /I n formed in gilt, and stored in temporary 
registers is used in gl5, thus saving operations. 
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Instruction 


Input to MRA 


In MRA Transfer 

from MRA 


gl£. 


ica tl 


* y n-l /:C n-l 


4y n-l /:t n-l 


irrrr 


yl+c 


* y n-l 




its 


tl 






ica 


t2 






imr 


yl+2+c 


* y n 


y 2 /l 


ica 


tl 


+y n-l y/I n-l 


2 / T 2 / T 
Wv/ti 


imr 


a5+c 


Xds n 


(y^/i +y 2 . /I , )ds 
n n n-1 n-1 n 






n-1 


n 


iad 


dl+12 


+ > (y 2 /I +y 2 _/I _ 

/ . J n / nVl' n-1 


)ds n )> (y n/ 1 n*Vl /1 n-l )dS n 






1 


1 


its dl+12 




n 

^ (y 2 / 1 +y 2 n A -i )ds 

/ . n n n-1 n-1 n 

1 


ict 


gll+2 


(Cycles back to gll+2) 


• 


isp 




(Address set by ita instruction at gll)« 



Due to symmetry of the web frame structure, 
dundant integrals are zero. Changes in slope due to 



the following re- 

V 



A9 




(+x) -yds 




Q -P 

(+x)-jds = q Q 




= 0 
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Vertical deflection due to m^ and p^. 




The sub-routines below perform the integration of moments caused 
by water pressure and weight loads e Note that two programs are necessary 
since h9 stations are used in the integration but only 25> x and y coordinates 
had been stored. The following program performs the integration of the 
moments along the starboard side to find the change in slope for one-half 
the web. 



A8 




m 

-f ds 




m _ _ 

n-l \ds 

I , } 2 
n-1 
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Instruction 

gl6, ica b3+c 
iad b3+2+c 
imr a5+c 
imr cl 

iad dl +50 



Input to MRA 



In MRA 



Transfer 
from MRA. 



+U1 /I 
n-1 7 n-1 



+m /I 
n n 



X ds 

A n 

X 1/2 



+m /I , 

n-1 n-1 

(m /I +m /I ) 

n n n-1 n-1 

(m /I +m /I )ds 

n 7 n n-1 7 n-1 n 

(m /I +m _/l )ds /2 

n n n-1 n-1 n 



n-1 n 

^ (m /I +m n /l _ )ds /2 y (m /I +m ../I _ )ds /2 

X .v n n n-1 7 n-1 n / X n n n-1 n-1 n 



its dl+50 



n 



A (m /I +m /I. )ds /2 
X . n n n-1 7 n-1 n 



The segment of the sub-routine that follows forms the integration 
of y m^/I to give the deflections in the x-direction around the starboard 
half of the web. 



Instruction 



Input to MRA 



In ~MRA Transfer 

from MRA 



gl7. 


ica yl+c 


+y n-l 


* y n-l 


imr 


b3+c 


X VaAi-l 


y n-l VXn-l 


its 


t3 




y n-l VlX-i 


ica 


yl+2+c 


+y n 


+y „ 


imr 


b3+2+c 


x m A 

n n 


y m / 1 
n n 7 n 


iad 


t3 


J n-1 n-1 7 n-1 


y m /I +y m ../I .. 

n n n n-1 n-1 7 n-1 


imr 


a5+c 


X ds 

^ n 


(y m /I +y _ m ,/l )ds 

n n 7 n J n-1 n-1 7 n-1 n 
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Instruction 



irar cl 



iad dl+52 



Input to MRA In MRA 

(continued) 



Transfer 
from MRA 



x 1/2 

n-1 



(y n m/ 1 * 7^1 » i/ 1 J ds /2 

n n n n-l n— _i_ n— 1 n 



+ / (y m /I +y n m n /l , )ds /2 
/ v J n n 7 n 17 n-1 n-1 7 n-1 n 7 



n 



(y m /I +y m /I )ds /2 
J n n 7 n J n-1 n-1 7 n-1 n 7 



its dl+52 



n 



(y m /I +y .m ./I .. )ds 

n n 7 n J n-x n-1 7 n-1 r 



To obtain the deflections in the y-direction on the starboard 
side due to applied loads, the integral of x m^/l was formed. 



Ay 



A 

/ xm 

t* 

0 



n 



x m _x _ mxds 
( n-1 n-1 . nn^ n 

' I + I ' 2 
n-1 n 



Instruction 

gl8, ica xl+c 
imr b3+c 
its t3 
ica xl+2+c 
irnr b3+2+c 



Input to MRA 


In MRA 


+x n 


+x 


n-1 


n-1 


ylH /I 

n-1 7 n-1 


x ,m ./I _ 
n-1 n-1 n-1 




X 




n- 


+x 


+x 


n 


n 


x m /I 


x m /I 


n in. 


n n n 


+ x ..m /I 
n-1 n-1 n-1 


(x ,m . n /l + x m /I ) 

n-1 n-1 7 n-1 n n 7 n 



Transfer 

MRA 



x ..m n /I .. 
n-1 n-1 7 n-1 



iad t3 
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Instruction 



Input to MRA 



In MRA 



Transfer 
from MRA 



imr a5+c 
imr cl 



(continued) 



* ds 

n 

x 1/2 



(x ..m n /I n +x m /I )ds 
n-1 n-1 n-1 n n n n 

(x _m _/l _+x m /I )ds (2 

n-1 n-1 n-1 n n n n 



n-1 



iad dl+5h + N )> (x n m /I +x m /i )ds /2 

/_ ^ n-1 n-1' n-1 n n n n' 






its dl+5U 



n 



/ (x _m ./I n +x m /I )ds /2 
Z » n-1 n-1 n-1 n n n' n 



ict gl6+2 



(Cycles back to gl6+2). 



To perform the integrations on the port side, the address changes 
if the x coordinates, y coordinates and girth segments ds must move up a 
column of registers,, This is accomplished by subtracting 2 from the address 
for each cycle performed. 

Formation of the change in slope due to applied loads (port side) 



AS 



O' 



n 



ds = 



25 



-7 m m _ ds 

(-E+ .”=!)— B- 

^ K I I /2 

n n-1 
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Instruction 



Input to MRA 



In MRA 



Transfer 
from MRA 



icr pbl 



(Sets cycle counter to number of points in a side) 



gl9. 


ica b3+pal + pal+c 


•ViX 


m ,/l , 
n-1' n-1 


iad 


b3+2 + pal+pal+c 


+ m /I 
n n 


(m /I +m /I ) 
n n n-1' n-1 


imr 


a5+pbl+pbl 


x ds 
* n 


(m /I +m n /l n )ds 
n n n-1' n-1 n 


imr 


cl 


xl/2 


(m /I +m n /l n )d.s /2 
n n n-1' n-1 n' 



its d 3+!?0 



25+n-l 

“7 



2£+n 

T 



iad d3+50 (m /I +m _/l _ )ds /2 \ n /I +m _/I _ )ds /2 

/ ' n n n-1 n-1 n' / n n n-1' n-1 n 



1 

25+n 



\ (m /I +m ,/l , )ds /2 

X . n n n-1 n-1 n 



Formation of horizontal deflection port side 

25+n 

Z 

25 



Ax 



0 25+n , 

/ m y tt 7 m _y , m y ds 

-r * - y 

a n-1 n 



Instruction 



Input to MRA 



g20, ica yl+pbl+pbl y n-1 
imr b3+pal+pal+c x m^ .^/i 
its tl 



In MRA 



n-1 



y _m _/I _ 

^n-l n-1' n-1 



Transfer 
from MRA 



y _m ./I _ 

n-1 n-1' n-1 



+y. 



n 



*7. 



n 



ica yl+pbl+pbl-2 
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Instruction 


Input to MRA 


In MRA 


Trans 

from MRA 




(continued) 






imr b3+pal+pal+2+c 


y m /I 
* n n 


y m /* 
n n n 




iad tl 


♦ViVA-i 


y m /I +y m ../I _ 

n n n ‘'n-l n -l n-1 


imr a^+pbl+pbl 


* ds n 


(y A *7 -i m 
n n n n-1 n- 


/I )ds 

•1' n-1 n 


imr cl 


Xl/2 


(y m /I +y n m 
n n n ‘'n-l n - 


/I )ds /2 
•l' n-1 n 



25>+n-l 

iad dl*S2 ^ 'Vn^Vn-A-l 14 / 



> (y m /I +y _m ./I , )ds /2 

/ „ n n n ^n-1 n-1' n-1 n' 

1 

its dl+52 



Formation of vertical deflection port side due to applied 



loads. 



Ay 



0 ' 2!?+n 

r ™ n X -7 ,VlVl . VnA 

-J -r ds - / , * — h- 

A n_1 n 



Instruction 


Input to MRA 


In MRA 


Transfer 
from MRA 


ica xl+pbl+pbl 


+x i 
n-1 


+X 1 

n-1 




imr b3+pbl+pbl+c 


+m . /I ~ 
n-1' n-1 


x n m _ /I . 
n-1 n-1 7 n-1 





,1 -./I -| 

n-l n-1 n-1 



its tl 



Instruction 



Input to MRA 



In MRA 



Transfer 
from MIA 



(continued) 



ica 


xl+pbl+pbl-2 


+x 

n 


+x 

n 


imr 


b3+pal+pal+2+c 


x m /I 
A rr n 


x m /I 
n n n 


iad 


tl 


VivA-! 


(y m /I +y _m /I _ ) 
n n n n-1 n-1 n-1 


imr 


a5+pbl+pbl+c 


* ds n 


(y rc /1+y , m /I )ds 

n n n n-1 n-1 n-1 n 


imr 


cl 


x 1/2 


AVVviv/W'V 2 




25>+n-l 






iad 


dl+5ii \ 

1 


(x m /I +x _m 
n n n n-1 n- 


.A-P'V 2 








2!?+n 








\ (x m /I +x _m _/l _ )ds /2 

/ v n n' n n-1 n-1 n-1 n' 



1 



its 



25+n 



(x m /I +x ..m 
n n n n-1 n- 



n 



The following program changes the addresses of the instructions 
in gl9 and g20 that utilize x, y and ds. The new addresses must move up 
a column of registers, 

OUT (Program goes into Whirlwind mode to change addresses). 
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Instruction 


Input to MRA 


In MRA 


Transfer 
from MRA 


ca 


gl9+2 


imr a5+pbl+pbl 


imr a£+pbl+pbl 




su 


c7 


-2 W.W. 


imr a£+pbl+pbl-2 




ts 


gl9+2 






imr a3+pbl+pbl-2 


ts 


g20+6 






imr a5>+pbl+nbl-2 


ts 


g20+17 






imr a3+pbl+pbl-2 


ca 


g20+3 


ica yl+pbl+pbl-2 


ica yl+pbl+pbl-2 




ts 


g20 






ica yl+pbl+pbl-2 


su 


c7 


-2 W.W. 


ica yl+pbl+pbl-U 




ts 


g20+3 






ica yl+pbl+pbl-U 


ca 


g20+lli 


ica xl+pbl+pbl-2 


ica xl+pbl+pbl-2 




ts 


g20+ll 






ica xl+pbl+pbl-2 


su 


C7 


-2 W.W. 


ica xl+pbl+pbl-Ij 




ts 


g20+ll|. 






ica xl+pbl+pbl-l; 


IN 


(Puts computer back in interpreted mode). 




ict 


gl9 








isp 











The above integrations provide the three equations for solving 
an unsupported ring frame of three degrees of indeterminancy. To solve the 
subject problem, the equations 



B 



f 



2 m(x-x )ds 
a 



B 



and 



li m(x-x^)ds 



B. 
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must be formed and solved simultaneously with the three previously formed. 
The designations and are the end points of the two stanchions, 

The sub-routine that follows determines the vertical deflection 
between B^ and. B^ of the starboard 



B 



B, 



Ay =. 



2 m (x-x ) 
n a 



'v m (x -x ) 



B, 



ds 



r m 1.x ■ 



m 



(x -x ) 



B, 



n-1 n-1 a \ ds 

1 i ^2 

n-1 



B, 



2 r 



Vl 

B. 



m (x -x ) . 
n n a ds 



n 



m i ( x ,-x )ds 
n-1 n-1 a n 



n-1 



Instruction 
ita g21+lfL 


Input to MRA 


In MRA Transfer 

from MRA 


g21, icr pa£ 


ica xl+pbl+c 


X 1 

n-1 


X 1 

n-1 


isu xl+pbh 


-x 

a 


(x ,-x ) 
n-l a 


its t3 




(x n-l" X } 
n-i a 


imr a5+pnli+c 


x ds 


( vr x a )ds 


imr cl 


x 1/2 


(x n-l' x a )ds/2 


its' tl 




( Vl’ x a )ds / 2 


Imr b3+pbli+c 


x (m/I) n-1 


(x -X )(m/l) ,ds /2 _ 

n-1 a ' n-1 n' n-1 




n-2 


n-1 






V 7 


iad dl+36 




/ , ( Vl- x a )(m/I) n-l dS n /2 




B i 


B 1 
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Instruction 



its dl+56 



Input to MRA In MRA 

(continued) 



Transfer 
from MRA 



n-1 



(x ..-x )(m/l) ds /2 
Z. x n-1 a n-1 n' 



B„ 



n 



ica xl+pbIi+2 +x 

isu xl+pblt -x 

its tU 

imr a5+pbli+c x ds 

imr cl X 1/2 

its t2 

imr b3+pbl+pbl+2+c x(m/l) 



n 



iad dl+56 



B +n-l 

—7 




+x 

n 

(x -x ) 
n a 



(x -x )ds 
nan 

(x -x )ds/2 
n a ' 



(x -x ) 
n a 



ds/2(x -x ) 
' n a 



(m/l) (x “x )ds/2 
n n a 



B^+n 
~ — 7 



(m/l) (x -x )ds/2 
' n n a ' 

B i B i 



(m/I) (x -x )ds/2 
x n n a 



V n 



its dl+56 



r \~ Z (m/I) (x -x )ds/2 
/ . n n a ' 

B, 



The following instructions perform the integration 
B 



n 

y(x-x ) 



0, 



ds = p r 



B, 





y (x -x ) 


£ 


n n a 


I 

n 



Is y n (x _-x ) ds 
n •'n-l n-1 a n 

2 I . 2 

n-1 
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Instruction 

ica tl 
imr d3+pbl|+c 
its tl 
imr yl+pbU+c 



Input to MRA 
(x n-l- X a )ds/2 

* VVi 



In MRA 



(x -x )ds /2 
n-1 a n 



Transfer 
from MRA 



xy 



n 



rl 



(x -x /I . )ds /2 
n-1 a n-1 n 

(x ,-x /I . )ds /2 
n-1 a! n-1 n / 

y (x -x )/l , ds /2 
n-1 n-1 a n-1 n 



B^+n-2 



B^+n-1 



iad dl +16 + N '> y . (x ,-x )/I ,ds /2 'N y n (x n -x )/l n ds /2 

X . n-1 n-1 a 7 n-1 n' / . n-1 n-1 a ' n-1 n' 

B, B 



B^+n-1 

“7 



its dl+16 
ica t2 

imr d3+pbIi+2+c 
its t2 



(x -x )ds /2 
n a n 



(1/D 



n 



B. 



(x -x )ds /2 
n a n 

(x-x )/l ds /2 
n a n n 



y -i (x ,-x )/I .ds /2 
X Xn-1 n ~l a n-1 n 



(x -x )/l ds /2 
n a n n 



imr yl+pbii+2+c x y 



n 



(x -X )y /i ds /2 
n a n n n 



B^n-1 



B^+n 



iad dl +16 t \ y . (x -x )/l ..ds /2 \ y (x -x )/l ds /2 

X ^ n- l n a n-1 n X Xn n a ' n n 

B. B. 



B_+n 




y (x -x )/I ds /2 
n n a n n 



its dl +16 
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The vertical deflection between the ends of the starboard stanchion 
due to the redundant force was found by 

B 

1 x(.x-x^, 

"-ds = (JqV 



B, 



/ 2 x(x-x ) —.x - (x _-x ) ds y (x -x ) ds 

a n-1 n-1 a n . J _n n a n 



V~> n-I n-1 
Za T n-1 



2 I 2 

n 



Instruction 



Input to MRA 



In MRA 



Transfer 



ica tl 


(x -x /I )ds /2 

n-1 a n-1 n 


(x -x )/l ds /2 
n-1 a n-1 n' 


imr xl+pblt+c 


* x n-l 


x n (x -x )/I _ds /2 
n-1 n-1 a ' n-1 n 




B^+n-2 


B + n -l 


iad dl +26 


' 7 V^n-Xa’XV 2 
B 1 


Zvi ( vr x a )/I „-iX 

B i 


its dl +26 




B +n-l 

1 

7 v/vpA 
B 1 


ica t2 


(x -x )/l ds /2 
n a n n 


(x -x )/I ds /2 
n a n n' 


imr xl+pbIt+2+c 


* X n 


x (x -x )/l ds /2 
n n a ' n u 



B +n-l 



B^+n 



iad dl +26 \ x , (x -x )/l _ds /2 \ x (x -x )/I ds /2 

X . n-1 n-1 a n-1 n' X n n a n . n 

B„ B, 



B^+n 



\ x (x -x )/I ds /2 
X . n n a n n 

B_ 



its dl +26 



13b 



The vertical deflection between the ends of the starboard stanchion 
due to the redundant column axial force R was found by 

B. 



x-x 



R J (x-x a ) ( j a )ds 
B„ 



Instruction 

ica t3 
imr tl 



(x -x ) 
n-1 a 



rV 


lx -x ) 
n a 


= L 


1 

n 


> MRA 




) 





n n-1 a n 

2 + I _ 2 

n-1 



In MRA 



(x -x ) 
n-1 a 



Transfer 
from MRA 



X ( Vr x a )/I n-l ds n /2 ( Vr x a ) 2 / I n-l ds/2 



B^+n-2 



B +n-l 



iad dl+36 + ^ (x -x ) 2 /l n ds /2 (x _-x ) 2 /I _ds /2 

Z n-1 a n-1 n Z ^ n-1 a ' n-1 n 

B_ B_ 



its dl+36 



B^+n 



(x -x ) 2 /I ds /2 
Z v n-1 a ' n-1 n 



B„ 



ica tU 
imr t2 



(x -x ) 
n a 

x (x -x )/l ds /2 
n a n n 



(x -x ) 
n a 

(x -x ) 2 /I ds /2 
n a ' n n 



B^+n-1 



B_ +n-l 



iad dl+36 -V" U -x ) 2 /I ds/. Y^(x -x ) 2 /I ds /2 

Z i. n ~- L a n-1 ^ z * n a n n 

B, B„ 



B^+n 



'S (x -x ) 2 /I ds /2 
/ „ n a n n' 



its dl+36 
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The vertical deflection between the two stanchion ends starboard 
side due to was performed below. 




'2 (x-x ) 

J ds = m 

B. 



(x -x ) ds (x _-x ) ds 
n a n n-1 a n 



n 



n-1 



Instruction 



Input to MRA 



In MRA 



Transfer 
from MRA 



ica tl 
iad t2 



(x _-x )/I ..ds /2 
n-1 a n-1 n 

4- (x -x )/I ds /2 
rr a n n 



(x -x )/l ..ds /2 

n-1 a n-1 n 

(x -x )/I _ds /2 

n-1 a ' n-1 n' 



+(x -x )/l ds /2 
n a n n' 



iad 



B 1+n -l 



dl+6 ■+ \ 7 (x --x )/I _ 

/ ^ n-1 a n-1 

B, 



+ (x -x )/l ds /2 
nan n 



B^+n 




( Vl- x a ,/l n-l dS n/2 



+(x -x )/l ds /2 
n a n n 



its dl+6 



ict g21+l 
isp 



B^+n 



B. 



[(x --x )/l - 

| n-1 a ' -n-1 



+ (x -x )/I "1 ds /2 
n a nl n 



(Cycles back through the stanchion integration 
program) . 



The following sub-routine performs the integrations around the 
port side which will give the total vertical deflection between the ends 
of the port stanchion. The entire sum should equal zero, since the stanchion 
is assumed incompressible. 
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The instructions that follow perform the integration 



B, 



B 



Ay 



B. 



m (x-x^) 



ds = 



In 






m 



^• ( V x b ) ( VrV 

n-1 



n 



Instruction 

ita g22+l8 
g22, icr pa6 
ica xl+pb5+c 

isu xl+pb5 

its t3 

imr a3+pb3+c 

imr cl 

its tl 

imr b3+pb6+pb6 



Input to MRA 



In MRA 



A n-1 

'*b 

x ds 
A n 

X 1/2 



* ^ rn / I ) n - 1 



B^+n-2 



Transfer 
from MRA 



- 

n-1 

Vi'^b 



(x n-l‘ x b ) 

( Vl" x b )dS n /2 

( Vl'*b )d V 2 

( m / I ) n -l (x n-l -x b )ds n /2 



B^+n-1 

7 



iad dl+58 (Vl)^ 1 <VlV < V 2 y %K-l" x b )/l n-l d V ' 2 

B„ B_ 



B +n-l 



its dl+58 



ica xl+pb^+2+c +x 



m (x .-x, )/I ,,ds /2 
n-n n-_ b n-1 n' 



B, 



n 



+x 



n 



■*b 



isu xl+pb5 



x n"*b 
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Instruction 


Input to MRA 


In MRA Transfer 




from MRA 




(continued) 




its tU 




V*b 


imr a5+pb5+c 


X ds 

n 


g 

CO 

> 


imr cl 


X 1/2 


XW 2 


its t2 




( VV ds n /2 


imr b3+pb6+pb6-2 


X (m/l) n 


^ m/l) n (x n*" x b )ds n /2 


B+n- 


•2 B +n 

x — 7 


iad dl+58 


(>"/i) n . 1 (v x b )ds / 2 / 


Wl) n U n -^)ie n /2 


B 3 


B 3 


V”. 


its dl+58 




/ ^nXW 2 

B 3 


OUT 






ca g22+l6 


imr b3+pb6+pb6-2 


imr b3+pb6+pb6-2 


ts g22+7 




imr b3+pb6+pb6-2 


su c7 




imr b3+pb6+pb6-li 


ts g22+l6 




imr b3+pb6+pb6-h 



IN 



The following section of the sub-routine performs the integral 



B. 



p d. 



y(x-Xh) 



ds = p, 



0 



V n-l ( Vl- x b ) . VvV 



B; 



n-1 



n 



ds 



n 



Instruction 



ica tl 
imr d3+pb5+c 
its tl 
imr yl+pb£+c 

iad dl+18 

its dl+18 
ica t2 

imr d3+pb5+2+c 
its t2 

imr yl+pb£+2+c 
iad dl +18 



Input to MRA 



In MRA 



1 rangier 
from MRA 



(x ,-x. )ds /2 
n-1 d n 






xy n-l 



( vr x b ,, V 2 

( Vl' x b )/I n-l d8 n ^ 

y n-l (x n-l‘ X b )/l n-l d V' 2 



B^+n-2 



B +n-l 

^ 7 



/. y n ( V X b ,/I n d3 n ' /2 3 r n-l ( Vl* X b )/I n-l ds n /2 

B„ B„ 



( V x b )ds n /2 



x 1/1 



n 



W. 



n 



B +n-l 



Zl ^XVXt/^n- 

B 3 



( VW 2 

( v\ )/I . d V ! 



(x “X, )/ . 
n o 



y (x -x, )/I ds /2 
n non n' 



B +n-l 



B f n 



X y _ (x _-x, )/l ,ds /2 x ; y (x -x, )/l ds /2 
X . J n-1 n-1 b n-1 n X . r. n V‘ n n 

B„ B- 



B^+n 



y (x -r )/l ds / 
X ' n n b n n' 



its dl +18 
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‘The following section of the sub-routine performs the integration 



B, 

nk 






x(x-x, 0 ) 



B 



ds = 



k 






B. 



X (x -x. ) , x _ (x ,-x, ) 
n n d ds n-1 n-1 b 



I 2 



n 



n-1 



Instruction 

ica tl 
imr xl+pb5+c 



Input to MRA 



In KRA 



Transit u 

from 



(x -x. )/I ds /2 (x -x. )/l ..ds /2 
n-1 b ' n-1 n n-1 o' n-1 n 



X x 



n-1 



x (x -x )/l ds /2 
n-1 n-1 d n-1 n 



B^+n-2 



B^+n-1 



iad dl+28 2l X n ( V x b )/I n ds/2 XL Vl'VfW/ 

B„ B. 



its dl +28 



B^+n-1 



B. 



x (x --X, )/I do ,/2 
n-1 n-1 b n-1 r 



ica t2 

imr xl+pb5+2+c 



W-W 2 



X x 



n 



(x -x )/I as /2 
no n n 

x (x -x, )/I ds /2 
n n d / n n 



B^+n-1 



B ? +n-l 



iad dl +28 / ^ x , (x -x )/l , ds /2 x , (x -x )/I ds /2 

/ . n-1 n b n-x n / „ n-1 n To " n it 

B. B. 



B^+n 



x (x -x )/l ds /2 
_s, n n o 7 n n 



its dl+28 



The following section of the sub-routine forms the vertical 



deflection caused by the axial stanchion force redundant S. 



B, 



Ay = S 



B. 



(x -x b )(x-x b ) 



II 



ds 



. ( Vl^b )2 1 ds n . 

°/ I 2 I , 2 

n n " 1 J 



Instruction 



Input to MRA 



In MRA 



Transfpu 
from MRA 



ica 


t 3 


Vi“*b 


1 

H 


imr 


tl 


X (x n-l -x b^ I n-l ds n / ^ 


(x ..-x,) 2 /! _ds /2 
n -1 b ' n -1 n' 






B +n -2 

V f - 


B +n-l 

t 0 


iad 


dl+l ;8 


/ (x n- 2 ~ x b ) ^n-S^n 72 


/ . (x n-l~ x b ) 7 l n-l ds n /2 






B 3 


B 3 








B +n-l 


its 


dl+li 8 




, < vi- x b ): 
B 3 


ica 


t h 


( VV 


(x n~ x b^ 


imr 


t 2 


x 'vV/VV 2 


^ x n -x b ^ 7 l n ds n 72 






B +n-l 

X 7 o 


B +n 

X r 9 


iad 


dl+U 8 


/ , ( Vr x b ) ‘ / 1 n-l dS n 72 


/ , ( Vl‘ X b ) 7 l n dS n /2 






B 3 


B 3 



n- 



B+n 



y (x -x ) 2 /I d s /2 
„ n b n n 



B, 



its dl+l |8 



The following section of the sub-routine finds the vertical 



deflection caused by the redundant moment m Q . 



B, 



by = m 



(x -V 



B, 



0 



ds = m„ 



B. 



' ( vV ds n . ( vrV 



n 



n-1 J 



ds 



Instruction 

ica tl 
iad t2 



iad dl+8 



its dl+8 



Input to MIA 
( Vr X b )/I n-l ds n /2 

K^ )/ 1 n d V 2 



B +n-l 

^ y 



B. 



In MRA 



Transfer 
from MRA 



( x n-l -x b > h„-l ^ s n / 2 

* W ^,/ 2 



B l" 



, (x n-l- x b )/I n-l 

* ( V x b )/I n] d V 2 



(x -x, )/I ds /2 
' n-1 b n-1 n 



B, 






B +n 

•3 



B, 



.<Vrt !/I rl 

4 ^-V 1 n]‘V 2 



ict g22+I 



(Cycles back to g22+l). 



The redundant stanchion forces affect the deflections for the 
web as a whole in much the same manner as the three basic redundant forces 
affect the deflection around the two stanchion sections. An analysis of 
Figure VI shows that the integrals of p^, q , and m for the traverse 

i-O* ^o» 0 

around the stanchion loop is the same as the integrals of R and S about 



the entire web girth. Therefore, the coefficients of the redundants which 
are the integrals may be used in several places . As an example the coeffi- 
cient of the deflection magnitude due to m^ between and is the same 
as the coefficient of the angular rotation magnitude due to the redundant 
R about the entire girth. Moment due to R is R(x-x ) and the rotation 

cl 

due to the moment is 




The deflection between B^ and B^ due to m^ is 



r 



nv 



(x ' x a ) 



• as. 



J 



The following instructions utilize the computed coefficients to 
fill similar values. 



Instruction 



Input to MRA 



In MRA 



Transfer 
from MRA. 



ica dl+6 



(x-x )/lds 
a 



B„ 



r 3 

+ I (x-x a ;/ld£ 



its dl+30 



B. 



J'' (x-x a )/ld 

B„ 



B 2 B 2 

y(x-x a )/ Ids + y(x-x a )/ Ids 



1 



B, 



ica dl+l6 



Instruction 



Input to MRA 



In MRA 



Transfer 
from MRA 



(continued) 



its dl +32 



ica dl +26 




its dl+3li 



ica dl+8 




(x-x, 0 )/lds 




y(x-x )/Ids 

3 . 




its dl+iiO 



ica dl +18 




its dl+li2 



B, 



J ° x(x-x a )/lds 

iL 




(x-x b )/Ids 




ica dl+28 



Instruction 



Input to MRA 



In MRA 



Transfer 
from MRA 



(continued) 



its dl+32 



ica dl+26 




its dl+3l| 



ica dl+8 




its dl+liO 



ica dl+l8 




its dl+U2 



J ^ x(x-x a )/Ids 





x(x-x )/lds 

cl 




(x-x^)/Ids 




ica dl+28 



lljli 



Instruction 



its dl+14; 



Input to HRA In MRA 

(continued) 



Trsnsfg-^ 
from MRA 




x(x-x cj )/lds 



isp 



The following sub-routine reduces the matrix of deflection 
equations by Crout's method and finds values for all redundants. The sub- 
routine, designated by S2, is a routine contained in the library of the 
computer laboratory. The library routines have been tested and tried 
and are known to be correct. For this reason only the method of entering 
the routine will be given. 



Instruction 



ita g23+£ 
OUT 

g23, SP S 2 
dl 

5 

1 



IN 



Input to MRA In MRA Transfer 

from MRA 



Address of first column. 

Order of matrix. 

Number of columns of constants 



isp 

The instructions below alter the vertical and horizontal forces 
to include the redundant forces p p , q^, R and S. Note that each redundant 
is applied as a load at its point of action. 



*2 ■ \ - »0 



V 2 ■ V 1 - % * E + S ' 



Instruction 


Input to MRA 


In MRA 


Transfer 
fror- MRA. 


ita g2li+25 










g2l>, ica e5 


+V 

n 




+v 

n 




isu dl+51| 


% 








its eh 


The 


following 


instructions form H-p. 

0 


P 

a' 

1 


ica e3 


+H 

w 




+H 

w 




isu dl+^2 


" P 0 




nj 

1 

O 




its e3 








H " p 0 
W u 




The 


following instructions form V-R 




ica e£+pbii 


+V 

n 




+v 

n 




isu dl+36 


-R 




V -R 
n 




its e3+pt»l| 








V -R 
n 


ica e5+pb5+pa5+pa5 


+V 

n 




+V 

n 




iad dl+56 


+R 




V +R 
n 




its e5+pbu+pa3+pa5 








V +R 
n 


ica e!?+pb6+pb6 


+V 

n 




+V 

n 




isu dl+5>8 


+S 




V +S 
n 




its e3+pb6+pb6 








V +S 

n 


ica ab6+pb6-pa6-pa6 


+V 




+V 




iad dl+58 


+s 




+V+S 




its pb6+pb6-pa6-pa6 








V+S 



The Whirlwind instructions that follow change the address of 



g2l+l6 and g2l|.+l8. 
Instruction 


Input to MRA 


In MRA 


Transfer 


OUT 

ca g2li+l6 


ica e£+pb6+pb6-2 


ica ep+pbo+pb6-2 


from MRA 


su c7 


-2 


ica e5+pb6+pb6-l| 




ts g2ii+l6 
td g2li+l8 
IN 

ict g2ii+l6 






ica e/+pb6+pb' 
e5+pb6+pb6-Ii 



The f oil owing sub-routine finds the stresses at each station due 
to the axial loads (p), the shear loads (q) and the bending moments (m). 
After the various stresses have been formed the direct stresses and the 
maximum compression and tension bending stresses are combined to form the 
total stresses at the flange and at the shell surface. The relationships 
listed below are utilized in the sub-program. 



p = H cos / - V sin / 

q = H sin / - V cos j i 

Direct stress 0 = p/A 

a 

Shear stress <7^ - q/A 
Bending Stresses = ny/ I 





m 



shell = m/I (d, T . +t ) 
in i» • • s 

flange^ - m/1 



Ib7 



The following instructions use the above relationships for the 
starboard side of the frame. 



Instruction 


Input to 


ita g27 




g2 S, ica c6 


+0.0 


its t8 




its t7 




icr pal 


n-1 
^ 


ica t8 




iad e5+c 


+v 


its t8 




imr al+c 


X sin / 


its tl 


n-1 


ica t7 


*z- 



In MRA 



Transfer 
from MRA 



+ 0.0 

+0.0 in t8 

+ 0.0 



n-1 




iad e3+c +H 

its t7 

imr a2+c X cos ^ 

-V sin / 




H 




H 



n 

"'y H cos / 
y * H cos <f> 




isu tl 



In MRA 



Instruction 



Incut to MIA 



(cent inusd) 



idv 


d2+c 


-r A 


p/A 


imr 


el 


X 8 


8 p/A 


its 


xl+c 







jVqm MR A. 



p/ A = cr 



n 



n 



ica t8 



imr a2+c 



its tl 



ics t? 



imr al 



x COS / 



n 

+ 'Z H 

X sin / 



+ > V 



II 

X 



V ccs $ 



n 



H 



n 



H sin <f> 



n 

/ V cos / 



n 



n 



n 



isu 


tl 


- ^ V cos <f> 


H sin / 


idv 


d2+c 


-7 A 


q/A 


imr 


el 


X 8 


8 q/A 


its 


sl+c 







q/A ~07 



Note that( 7 ^ was formed using the entire area as effective in 
shears. This is an error which should be corrected since on].. 7 the area 
of the web is effective in shear. 



-d 



II. A. 



“ d N.A. 



ics dl+c 



Ih9 



Instruction 



Input to MRA In MRA Transfer 

from HRA 



(continued) 



isu 


b5>+c 


-t 

s 


-'■wv 


imr 


Cl5 


X 12 


- 12(d N,A.*V 


imr 


b3+c 


X m/I 


-12 n /l(d N-A _*t s ) 


imr 


el 


X 8 


-96 V«d N>A _-t s ) - - CQ 


iad 


xl+e 


+ 

a 




its 


a5+c 




<7T 


ica 


bl+c 


+d 

w 


+d 

w 


iad 


bii+c 


+t f 


Vh 


isu 


dl+c 


~Va. 


^ + 
1 

■*— H 

• 

> 

• 


imr 


cl5 


X 12 




imr 


b3+c 


X m/I 


12 ■nA(V t f' d H.A. ) 


imr 


el 


CO 


96 r»/I(d v * V d N>A> ) . <£ 


iad 


xl+c 


+ 

a 




its 


e5+c 




a~ f 


ict 


g2$+3k 







The following instructions calculate the stresses around the 



port side 



i5o 



Instruction 


Input to MRA 


In MRA 


Transfer 
from MRA 


icr pbl 


n-l 


n-l 




g26, ica t8 












n 




iad e5+pal+pal+c 


+v 


£. v 


n 


its t8 




n 


Zv 


imr al+pb2+pb2 


X sin / 


sin V 


n 


its tl 


n-l 


n-l 


+sin $ / V 


ica t7 


+ !> , H 










n 




isn e3+pb2+pb2 


+H 




n 

^ » 


its t7 




n 




iror a2+pb2+pb2 


X cos (f) 


cos <{>/ H 




iad tl 


n 

+sin V 


n 

cos / ^ H+sin 


P-i 


idv d2+pb2+pb2 


-f- A 


p/A = <7^ 




its xl+pal+pal+c 






p/A —o~~ 



a 



Instruction 


Input to MRA 


In MRA 


Transfer 
from MRA. 




(continued) 








n 


n 




ica t8 




Zv 








n 




imr a2+pb2+pb2 


X cos / 


cos / V 


n 

\ 


its tl 


n 


n 


cos / / 


ica t7 












n 

* 




imr al+pb2+pb2 


X Sin / 


sin / / , H 






n 


n 

— 7 


n 


isu tl 


-cos / y v 


+(sin / / ^ H-cos / / 


idv d2+pb2+pb2 


-j-A 


> 

I! 




its sl+pal+pal+c 








ics dl+pb2+pb2 


*" d N.A. 


" d N.A. 




isu b£+pb2+pb2 


-t 

s 


~ (d H.A. +t s > 




imr b3+pal+pal+c 


X m/I 


- n / l(d N.A. + V 


■ 


isu xl+pal+pal+c 


-<r 

a 


-( + ) = C7 ' 

a sin s 








its a5+pal+pal+c 
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Instruction Input to MRA In MRA Transfer 

from MRA. 

(continued) 



ics dl+pb2+pb2 


“ d N.A„ 


_d H.A. 


iad blj+pb2+pb2 


+t f 


' d N.A. +t f 


iad bl+pb2+pb2 


+d 

w 


( W <J K.A. ) 


imr b3+pal+pal+c 


X m/I 


= frii 


iad xl+pal+pal+c 






its e5+pal+pal+c 




<?T 


OUT (The 

the 


following instructions 
above program). 


change the addresses of 


ca g26+3 


imr al+pb2+pb2 


imr al+pb2+pb2 


o- 

o 

2 

CO 


-2 


imr al+pb2+pb2-2 


ts g 26+3 




imr al+pb2+pb2-2 


ts g26+17 






ca g26+7 


iad e3+pb2+pb2 


iad e3+pb2+pb2 


su c7 


-2 


iad e3+pb2+pb2-2 


ts g 26+7 




iad e3+pb2+pb2-2 


ca g26+9 


imr a2+pb2+pb2 


imr a2+pb2+pb2 


su c7 


-2 


imr a2+pb2+pb2-2 


ts g 26+9 




imr a2+pb2+pb?-2 


ts g26+lU 




imr a2+pb2+po2-2 


ca g 26 +ll 


idv d2+pb2+pb2 


idv d2 +pb2+pb2 


su c7 


-2 


idv d2+pb2+pb2-2 



ts g 26 +ll 



idv d2+pb2+pb2-2 
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Instruction 


Input to MRA 


In MRA 


Transfer 










from MRA 






(continued) 






ts 


g26+l8 






idv d2+pb2+pb2-2 


ca 


g26+20 


ics dl+pb2+pb2 


ics dl+pb2+pb2 




S11 


c7 


-2 


ics dl+pb2+pb2-2 




ts 


g26+20 






ics dl+pb2+pb2-2 


ts g26+25 






ics dl+pb2+pb2-2 


ca 


g26+21 


isu b5+pb2+pb2 


isu b5+po2+pb2 




su 


c7 


-2 


isu b5+pb2+pb2-2 




ts 


g26+21 






isu b5+pb2+pb2-2 


ca 


g26+26 


iad bit+pb2+pb2 


iad bli.+pb2+pb2 




su 


c7 


-2 


iad bIi+pb2+pb2-2 




ts 


g26+26 






iad bli+pb2+pb2-2 


ca 


g26+27 


iad bl+pb2+pb2 


iad bl+pb2+pb2 




su 


c7 


-2 


iad bl+pb2+pb2-2 





ts g26+27 
IN 

ict g26 
g27, isp 



iad bl+pb2+pb2-2 
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